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SUMMARY 


Experiments  were  conducted  with  two  Delta  wings  and  a  streamlined 
body  with  lifting  surfaces  to  Investigate  the  evolution  of  turbulent 
trailing  vortices  In  stratified  and  unstratified  water.  The  vortex 
trajectories  were  determined  as  a  function  of  the  relative  depth  D/bQ, 
normalized  time  VQt/b0,  and  the  stratification  parameter  Nb0/VQ.  The 
results  have  shown  that  the  vortices  rise  only  to  a  finite  height  as 
they  gradually  demise  under  the  Influence  of  turbulence,  sinusoidal 
instability,  and  the  vortex  breakdown.  The  effect  of  stratification 
Is  to  reduce  the  lifespan  of  vortices  and  the  maximum  height  attained 
by  them.  Various  approximate  analyses  have  been  examined  and  new  ones 
have  been  proposed  to  determine  the  lifespan  of  the  vortices.  Finally, 
an  exploratory  numerical  analysis  of  the  two-dimensional,  unsteady, 
laminar  vortices  has  been  carried  out  with  encouraging  results. 
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I.  INTRODUCTION 


Trailing  vortices  formed  by  a  lifting  surface  traveling  through  a 
stratified  or  unstratified  fluid  medium  have  been  the  subject  of  recent 
interest  primarily  because  of  their  importance  in  the  determination  of 
safe  aircraft  separation  distances  for  terminal  areas.  In  recent  years 
this  interest  has  extended  to  the  trailing  vortices  generated  by  the 
control  planes  of  submerged  bodies  such  as  a  submarine  partly  because  of 
the  effect  of  these  vortices  on  the  motion  of  the  submerged  body  itself 
and  partly  because  of  the  interaction  of  these  vortices  with  the  surround¬ 
ing  fluid  medium. 

The  trailing  vortices  may  be  generated  by  any  number  of  lifting  sur¬ 
faces  such  as  an  airplane  wing,  a  missile,  single  or  double  Delta  wings, 
or  the  control -surfaces  of  a  submerged  body  (sail  planes,  rudder,  and 
stem).  The  initial  organization  or  roll-up  of  a  trailing  vortex  from 
the  vorticity  generated  by  the  lifting  surface  is  very  rapid  and  is 
generally  complete  within  the  space  of  several  lifting- surface  spans. 

Whereas  a  single  vortex  is  extremely  stable,  one  or  more  vortices  in 
close  proximity  exhibit  a  number  of  complex  interactions  and  give  rise  to 
new  phenomena  which  tend  to  accelerate  the  demise  of  the  vortices  in¬ 
volved.  When  a  lifting  surface  starts  from  rest,  a  starting-vortex  or 
a  Horseshoe  vortex  is  formed  since  a  vortex  must  be  connected  either  to 
itself  or  to  a  rigid  boundary.  A  Horseshoe  vortex  is  also  formed  when  a 
lifting  surface  accelerates  or  decelerates  rapidly,  or  undergoes  a  radical 
change  in  the  angle  of  attack.  When  a  control  surface  (e.g.,  a  Delta 
wing)  starts  from  rest  and  moves  with  a  constant  velocity,  the  Horseshoe 
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vortex  closes  on  itself,  breaks  away  from  the  rest  of  the  trailing  vor¬ 
tex,  links  rapidly,  and  forms  a  vortex  ring.  While  the  first  vortex 
ring  undergoes  highly  complex  deformations,  a  second  vortex  ring  is 
formed  from  the  Horseshoe  vortex  connecting  the  tail  end  of  the  trailing 
vortex.  This  process  does  not  continue  forever  for  a  number  of  reasons. 
The  speed  of  the  generating  surface  is  much  greater  than  the  speed  of 
advance  of  the  vortex  rings.  Thus,  the  length  of  the  vortex  trail  in¬ 
creases  rapidly  and  the  motion  of  the  vortices  becomes  dominated  by  the 
mutual  interaction  of  the  vortices  rather  than  by  the  horseshoe  vortices 
at  their  tail  end.  The  motion  of  each  new  vortex  ring  differs  from  the 
one  shed  previously  partly  because  of  the  continued  aging  or  diffusion 
of  vorticity  (i.e.,  a  weaker  ring  is  formed  at  succeeding  steps)  and 
partly  because  the  motion  of  each  new  ring  is  controlled  by  all  the  rings 
shed  previously  and  by  the  ever-increasing  vortex  trail.  These  phenomena 
lead  to  the  rapid  demise  of  vortex  rings.  In  fact,  shortly  after  the 
start  of  the  motion  of  the  lifting  surface  the  horseshoe  vortex  and  the 
vortex  rings  play  a  very  minor  role  on  the  evolution  of  the  trailing 
vortex. 

The  trailing  vortices  are  subjected  to  three  known  decay  mechanisms: 
(i)  Aging  of  vortices  due  to  eddy  viscosity;  (ii)  sinusoidal  instability 
leading  to  the  linking  of  vortices  (known  as  the  Crow  instabil ity) ;  and 
(iii)  vortex  breakdown. 

These  mechanisms  can  act  independently  or  in  combi  nation.  The  viscous 
and  turbulent  diffusion  of  vorticity  lead  to  the  gradual  dissipation  of 
the  trailing  vortices.  Clearly,  this  phenomenon  is  strongly  influenced 
by  the  initial  vorticity  distribution  and  the  ambient  turbulence.  Further 
more,  the  proximity  of  other  vortices,  free  surface,  rigid  surfaces 
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(e.g.,  the  submerged  body,  ground,  etc.),  ambient  shear,  stratification 
of  the  fluid  medium,  etc.,  affect  the  demise  of  the  vortices.  Clearly, 
the  aging  of  the  vortices  is  ever  present  and  affects  the  evolution  of  the 
trailing  vortices  regardless  of  whether  the  Crow  instability  and/or 
vortex  bursting  take  place  or  not. 

Breaking  up  by  linking  or  sinusoidal  instability  involves  the  growth 
of  a  rather  regular  sinusoidal  pattern,  usually  symmetrical  for  the  trail¬ 
ing  vortex  pair.  This  is  commonly  referred  to  as  Crow  linking.  The 
time  to  link  (the  formation  of  a  series  of  vortex  rings  and  Horseshoe 
vortices)  has  been  found  to  depend,  statistically,  on  atmospheric  turbu¬ 
lence  within  the  inertial  subrange  of  eddy  sizes  (a  size  range  from  a  few 
cm  to  some  hundreds  of  meters).  The  scale  of  linking  varies  around  nine 
times  the  vortex  spacing.  Thus,  it  is  the  eddy  sizes  in  the  range  of  one 
to  ten  times  the  span  or  vortex  spacing  which  has  the  main  effect  on  link¬ 
ing.  The  formation  of  vortex  rings  by  linking  (distinctly  different  from 
those  generated  at  the  start  of  motion)  lead  to  the  rapid  demise  of  the 
vortex  trail.  Consequently,  one  would  want  to  accelerate  the  formation 
of  links  (e.g.,  by  periodically  changing  the  angle  of  attack  of  the  control 
surfaces)  if  one  is  interested  in  decreasing  the  hazard  posed  by  the  trail¬ 
ing  vortices  to  other  aircraft  or  in  decreasing  the  energy  transferred 
to  the  fluid  medium. 

The  vortex  breakdown  is  essentially  a  single- vortex  phenomenon,  in  the 
sense  that  it  is  not  necessarily  connected  with  the  induced  flows  asso¬ 
ciated  with  the  other  vortex  of  the  pair.  There  are,  however,  frequent 
occurrences  of  core  bursting  in  conjunction  with  the  sinuous  instability, 
in  which  case  the  sinuous  vortex  deformation  appears  to  induce  core  burst¬ 
ing.  The  effect  of  vortex  breakdown  is  to  flatten  the  velocity  profile  by 
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significantly  reducing  the  peak  velocity  (or  the  velocity  of  numerous 
peaks  of  several  smaller  vortices)  under  the  center  of  the  vortex,  but 
without  major  effect  on  the  speeds  away  from  the  vortex  center.  A 
definite  vortex  flow  remains  in  all  cases,  indicating  that  core  bursting 
does  not  totally  eliminate  the  vortex.  In  fact,  the  persistence  of  the 
outer  flow  circulation  suggests  that  vortex  breakdown  may  not  signifi¬ 
cantly  alleviate  the  vortex  hazard  and  the  energy  transfer  to  the 
surrounding  medium.  The  understanding  of  the  vortex-breakdown  and  the 
determination  of  how  much  circulation  remains  after  breakdown  remain  as 
critical  and  vexing  questions. 

The  trailing  vortices  move  downward  (downwash)  when  generated  by  a 
lifting  surface  with  a  positive  angle  of  attack,  as  in  the  case  of  an 
airplane,  or  upward  (upwash)  when  generated  by  a  sail  plane  with  a  nega¬ 
tive  angle  of  attack,  as  in  the  case  of  a  submarine.  In  either  case,  the 
motion  of  the  vortices  is  determined  by  the  mutual  induction,  vorticity 
distribution,  the  three  types  of  the  instabilities  or  demise  mechanisms 
cited  above,  stratification  effects  of  the  fluid  medium,  ambient  flow 
turbulence,  shear  and  currents  present  in  the  fluid  medium,  and  the 
proximity  effects  (e.g.,  free  surface,  bottom,  ground,  other  bodies,  etc.). 
An  inviscid  flow  analysis  based  on  concentrated  line  vortices  can  help 
to  predict  the  early  stages  of  motion  of  the  vortices.  However,  the 
later  stages  of  the  motion  become  increasingly  dominated  by  the  factors 
cited  above  and  it  is  not  possible  to  examine  the  evolution  of  the  trail¬ 
ing  vortices  without  taking  into  consideration  the  effect  of  viscosity, 
turbulence,  stratification,  mutual  induction,  proximity,  and  all  the  other 
phenomena  which  arise  as  a  consequence  of  these  fundamental  variables. 
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It  must  also  be  pointed  out  that  whereas  the  trailing  vortices  generated 
by  an  aircraft  are  under  each  other's  influence  (in  addition  to  the 
external  non-vortical  influences),  the  trailing  vortices  generated  by  a 
large  submerged  body  may  also  be  under  the  influence  of  a  complex  system 
of  image  vortices  imbedded  in  the  body. 

Density  stratification  within  the  fluid  medium  is  of  great  importance 
in  the  migration  and  ultimate  demise  of  trailing  vortices.  The  atmos¬ 
phere  is  less  dense  at  higher  altitude  while  the  oceans  exhibit  various 
patterns  of  stratification  due  to  local  temperature,  salinity,  and  depth. 
Generally,  a  vortex  migrates  with  the  fluid  in  which  it  was  born,  retain¬ 
ing  most  of  the  original  fluid  in  a  so-called  recirculation  cell.  The 
upwash  of  the  vortices  in  the  ocean  environment  into  a  layer  of  lesser 
density  gives  rise  to  a  buoyancy  which  is  in  the  opposite  direction  to 
the  velocity  resulting  from  the  mutual  induction.  Furthermore,  the  motion 
of  the  recirculation  cell  produces  oppositely-signed  vorticity  in  the 
surrounding  medium.  This,  in  turn,  gives  rise  to  the  entrainment  and 
detrainment  of  fluid  from  the  recirculation  cell.  Altogether,  these 
additional  consequences  of  stratification  help  to  reduce  not  only  the  life 
of  vortices  but  also  the  height  to  which  the  vortices  rise  asymptotically. 
One  may  also  surmise  that  the  stratification  will  tend  to  change  the 
occurrence  and  the  resulting  characteristics  of  eddy  diffusion.  Crow 
instability,  and  core  bursting.  It  is  evident  from  the  foregoing  that 
the  phenomena  associated  with  the  generation,  migration,  interaction, 
and  demise  of  trailing  vortices  are  very  complex.  As  the  review  of  the 
previous  investigations  will  reveal,  it  is  not  yet  possible  to  devise 
mathematical  models  with  which  some  or  all  of  the  characteristics  of  the 
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phenomena  can  be  predicted.  As  in  many  other  flow  situations,  the  role 
played  by  turbulence  and  the  quantification  of  the  turbulence  stresses 
pose  difficult  questions.  The  gradual  understanding  of  the  behavior  of 
trailing  vortices  must  rely  for  the  time  being  on  carefully  conducted 
experiments  without  losing  sight  of  the  fact  that  one  must  ultimately 
devise  a  model  to  acquire  a  power  of  prediction  as  well  as  a  means  to 
reduce  the  number  of  experiments  for  model  verification. 

It  is  clear  from  the  foregoing  that  the  reason  for  developing  as  much 
an  understanding  as  possible  of  the  behavior  of  trailing  vortices  Is  to 
permit  design  and  improvement  of  an  operational  system  which  properly 
balances  safety  and  efficiency  in  the  motion  of  vortex-generating  bodies 
and  also  to  aid  in  the  development  of  devices  to  enhance  vortex  decay, 
limit  vortex  upwash,  and  minimize  energy  transfer  to  the  surrounding 
medium.  The  effort  towards  the  understanding  of  some  of  the  underlying 
phenomena  forms  the  basis  of  the  present  investigation. 
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II.  REVIEW  OF  PREVIOUS  INVESTIGATIONS 


Vortices  and  vortex  wakes  have  become  a  major  theme  of  aerodynamics 
research  since  the  advent  of  the  large  aircraft.  The  objectives  have  been 
to  predict  the  trajectories  and  lifetimes  of  trailing  vortices,  in  order 
to  relieve  the  constraint  on  flight  operations  imposed  by  the  vortex 
hazard.  Much  of  the  progress  made  during  the  past  two  decades  was  dis- 
cussed  at  the  Symposium  on  Aircraft  Wake  Turbulence  and  Its  Detection 
[Ref.  1],  held  in  Seattle  in  1970,  and  at  the  Aircraft  Wake  Vortices 
Conference  [Ref.  2],  held  in  Cambridge,  MA  in  1977. 

The  fundamental  questions  regarding  the  evolution  of  the  trailing 
vortices  may  be  summarized  as  follows:  (i)  how  does  the  vorticity  gener¬ 
ated  by  the  lifting  surfaces  organize  itself  into  two  or  more  centers  of 
vorticity  behind  and  to  each  side  of  the  lifting  surface?  (ii)  how  do 
the  discrete  centers  of  vorticity  interact  and  affect  the  stability  of 
each  other?  (iii)  how  does  a  vortex  demise  and  how  does  the  demise  of  a 
single  vortex  differ  from  that  of  a  vortex  pair?  And  finally,  (iv)  what 
is  the  effect  of  stable  stratification  on  the  generation,  roll-up,  inter¬ 
action,  stability,  and  demise  of  vortices? 

A.  ROLL-UP  OF  VORTICITY 

It  is  a  well  known  fact  that  a  vortex  sheet  is  trailed  from  a  finite- 
aspect-  ratio  lifting  surface  as  a  consequence  of  the  nonuniform  spanwise 
lift  on  the  wing.  The  roll-up  of  this  sheet  into  discrete  vortices  as 
a  result  of  a  convective  motion  was  recognized  as  early  as  1907  by 
Lanchester  [Ref.  3].  Since  then  there  have  been  numerous  studies  to 
describe  this  complicated  phenomenon.  The  first  mathematical  model  to 


18 


describe  the  Invlscld  structure  of  the  vortex  wake  was  proposed  by  Prandtl 
[Ref.  4].  He  based  his  model  on  the  conservation  of  mechanical  energy 
and  found  that  the  vortex  tube  radius  for  an  elliptically  loaded  wing 
is  approximately  8  percent  of  the  wing  span.  In  this  analysis  the  choice 
of  the  swirl  velocity  distribution  is  somewhat  arbitrary,  being  only 
constrained  by  circulation  at  large  radius  from  the  vortex  center  and 
the  integral  constraint  on  kinetic  energy.  In  other  words,  the  model 
does  not  address  itself  to  the  determination  of  the  final  structure  of 
the  wake  for  a  given  initial  vorticity  distribution. 

The  second  model,  which  overcomes  some  of  the  arbitrariness  inherent 
in  the  Prandtl  model,  was  developed  by  Betz  [Ref.  5].  Betz’s  model  is 
based  on  four  integral  invariants  of  an  incompressible  two-dimensional 
bounded  vorticity  distribution.  The  Betz  model  has  now  been  checked 
against  measurements  made  behind  model  and  full-scale  aircraft  and  it 
*  has  been  found  that  the  model  predicts  distributions  which  are  in  far 
better  agreement  than  the  Prandtl  model  [Ref.  4]. 

In  more  recent  times  the  advent  of  numerical  methods  and  computer 
gave  impetus  to  the  use  of  discrete  vortex  models.  In  most  numerical 
analyses,  the  attached  boundary  layers  on  the  wing  surfaces  are  repre¬ 
sented  by  bound-vortex  sheets  while  the  separated  free-shear  layers  are 
modeled  by  free-vortex  sheets.  The  free-vortex  sheets  join  the  bound- 
vortex  sheets  along  the  separation  lines  which  are  known  a  priori  for 
wings  with  sharp  edges.  Furthermore,  it  is  assumed  that  the  vortex 
breakdown  (see  e.g..  Refs.  6-8)  occurs  far  downstream  so  that  the 
primary  core  size  and  its  variation  in  the  vicinity  of  the  wing  are 
negligible.  This  assumption  limits  the  angle  of  attack  (dependent  on 
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the  aspect  ratio  of  the  wing)  at  which  the  Invlscld  model  Is  applicable. 
Furthermore,  the  flow  outside  the  bound-vortex  sheet(s)  and  the  free- 
vortex  sheet(s)  Is  assuned  irrotatlonal .  Within  these  limitations,  the 
resulting  potential  flow  model  represents  the  main  features  of  the  real 
flow  to  a  high  degree  of  accuracy.  The  literature  contains  several  steady 
and  unsteady  inviscid-flow  models  with  various  degrees  of  limitations 
and  drawbacks. 

The  first  group  of  models  uses  s lender-body  theory  and  conical  flow 
assumption  [Refs.  9-12].  These  models  satisfactorily  predict  the  pressure 
distribution  over  the  front  portion  of  the  wing  surface.  In  the  rear 
portion,  the  models  fail  to  predict  a  satisfactory  pressure  distribution 
because  the  Kutta-Joukovsky  condition  is  violated  at  the  trailing  edge. 
Such  models  are  limited  to  simple  delta  planforms.  The  second  group  of 
models  uses  a  nonlinear  discrete- vortex  method  [Refs.  13-20].  In  these 
models,  the  bound  vortex  sheet  and  the  free- vortex  sheets  are  approximated 
by  a  set  of  concentrated  vortex  lines.  The  bound  vortex  sheet  is  replaced 
by  a  bound-vortex  lattice,  while  the  free-vortex  sheet  is  replaced  by 
segmented  free-vortex  lines  (in  the  case  of  steady  flow)  or  by  growing 
free-vortex  lattice  (in  the  case  of  unsteady  flow).  The  boundary  condi¬ 
tions  are  satisfied  at  certain  control  points  on  the  bound-  and  free- 
vortex  system  using  an  iterative  technique. 

Although  the  discrete  vortex  model  has  worked  well  for  many  years  for 
attached  flow  problems,  when  vortex- type  separation  from  leading  edges 
and/or  tips  occurs,  the  free-vortex  system  lies  close  to  the  lifting- 
surface  bound  vortices  and  the  results  are  found  to  be  sensitive  to 
variations  in  the  shapes  of  the  vortex  elements  and  the  relative  lengths 
of  the  vortex  segments  (see  e.g..  Ref.  21). 


The  third  group  of  models  employs  doublet-panels  (see  e.g..  Refs. 
22-23).  In  this  method,  the  wing  and  Its  free-vortex  sheets  are  divided 
into  networks  of  quadrilateral  panels.  Each  panel  of  the  networks 
representing  the  wing  has  a  biquadratic  local  doublet  distribution  and 
a  bilinear  local  source  distribution.  Source  and  doublet  splines  are 
used  to  express  the  distributions  of  singularities  on  the  networks  in  terms 
of  discrete  values  of  singularity  strength  at  certain  standard  points  on 
each  network.  The  boundary  conditions  and  continuity  of  singularity 
strengths  across  abutting  networks  are  enforced  at  certain  standard  points 
on  each  network.  The  results  of  this  method  are  generally  good  when  the 
solution  converges  (see  e.g..  Refs.  24-25).  Apparently,  the  difficulty 
in  obtaining  convergence  is  due  to  the  failure  in  satisfying  the  con¬ 
tinuity  of  the  derivatives  of  the  doublet  strength  across  abutting  net¬ 
works.  This  is  equivalent  to  the  existence  of  concentrated  vortex  lines 
between  abutting  networks.  Additional  related  works  are  discussed  in 
Ref.  26  in  connection  with  the  AGARD  Symposium  on  High  Angle  of  Attack 
Aerodynami cs . 

The  fourth  group  of  models  employs  a  nonlinear  hybrid  vortex  method 
(see  e.g..  Refs.  27-28).  In  this  method,  conti nuous-vorti city  and  vortex¬ 
line  representations  of  the  wing  and  its  separated  free-shear  layers  are 
used.  Continuous  vorticity  is  used  in  the  near-field  calculations  while 
discrete  vortex  lines  are  used  in  the  far-field  calculations.  The  calcu¬ 
lated  net  surface  pressure  and  section  normal-force  coefficients  for 
wings  with  tip  separation  are  in  good  agreement  with  the  experimental 
data  and  the  nonlinear  discrete-vortex  method. 

The  roll-up  calculations  do  not  give  details  of  the  structure  of  the 
tip  vortex  to  be  used  as  initial  conditions  in  a  viscous  solution.  Often 
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the  approximate  solution  of  Kaden  [Ref.  29]  for  the  unsteady  roll-up  of 
an  infinite  sheet  is  used  as  the  basis  for  the  prediction  of  the  decay 
of  laminar  trailing  vortices  [Ref.  30].  The  most  obvious  drawback  of 
the  existing  models  is  the  lack  of  a  realistic  model  of  the  leading- 
edge  core  and  its  feeding  sheet.  Modeling  of  the  leading-edge  core  and 
its  feeding  sheet  had  been  first  introduced  by  Brown  and  Michael  [Ref. 
31].  This  model  has  been  improved  by  Mangier  and  Smith  [Ref.  32]  and 
Smith  [Ref.  33].  Additional  improvements  have  been  made  by  Moore  [Ref. 
34]  who  devised  a  method  to  collapse  the  tightly-wound  inner  portion  of 
the  vortex  sheet  into  a  single  tip  vortex.  This  was  done  primarily  to 
avoid  the  numerical  difficulties  associated  with  the  crossing  of  the 
vortex  sheets  and  with  the  spurious  interaction  between  the  successive 
turns  of  the  inner  portions  of  the  vortex  sheet.  Clearly,  this  sacri¬ 
fices  the  detail  of  the  tip  vortex  structure.  The  latest  efforts  regard¬ 
ing  the  calculation  of  the  roll-up  process  use  the  method  of  rediscreti¬ 
zation  of  shear  layers.  This  method  has  been  devised  by  Fink  and  Soh 
[Ref.  35]  and  extensively  used  by  Sarpkaya  and  Shoaff  [Ref.  36].  The 
numerical  models  based  on  the  inviscid  flow  theory  do  not  account  for 
the  diffusive  effects  of  viscosity  and  turbulence.  Thus,  they  are  not 
capable  of  predicting  the  evolution  of  the  trailing  vortices. 

B.  TRAILING  VORTICES 

The  vortex  sheets  roll-up  asymptotically  and  it  may  be  several  wing 
spans  behind  the  aircraft  before  the  fully  rolled-up  condition  is 
reached.  With  the  high  Reynolds  numbers  encountered  on  lifting  surfaces 
in  free  flight,  it  is  likely  that  their  boundary  layers,  and  thus  their 
trailing  vortex  sheets,  are  turbulent.  At  finite  Reynolds  numbers. 
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vorticity  diffuses,  and  the  sheet  becomes  a  vortex  layer  of  finite,  non¬ 
zero  thickness.  Successive  turns  in  the  centermost  portions  of  the 
rolled-up  layer  merge,  the  spiral  structure  disappears  and,  if  the  radius 
of  the  vortex  core  is  greater  than  the  diffusive  length  scale,  an  equi¬ 
librium  structure  ensues.  Thus,  one  anticipates  a  multi-structured 
core,  the  inne;most  part  of  which  has  effectively  reached  a  state  of 
equilibrium  and  surrounding  this  an  annulus  in  which  the  turns  have  merged 
or  are  in  various  stages  of  merging,  and  an  outer  region,  beyond  which 
discrete  turns  of  the  spiral  remain  evident. 

Dynamics  of  the  flow  require  axial  velocities  to  be  present  in  trailing 
vortices  [Ref.  37]  and  these  can  have  an  important  effect  on  the  stability 
of  the  trailing  vortex  systems  and  play  an  important  role  in  the  occurrence 
of  vortex  breakdown.  It  is  also  a  well-known  fact  that  the  occurrence  of 
vortex  breakdown  in  trailing  vortices  may  be  delayed  or  enhanced  by  the 
acceleration  or  deceleration  of  the  ambient  flow  due  to  the  proximity  of 
streamlined  axisymmetric  bodies  whose  control  surfaces  may  have  given  rise 
to  the  development  of  trailing  vortices.  Thus,  it  is  important  to  make 
a  distinction  between  the  trailing  vortices  of  an  aircraft  and  those  of  a 
submerged  body.  In  the  case  of  the  former,  the  trailing  vortices  are  not 
affected  by  the  proximity  of  the  aircraft  body  partly  because  the  length 
of  the  aircraft  in  terms  of  the  wing  span  is  relatively  small  and  partly 
because  the  wing  span  is  many  times  that  of  the  body  diameter.  In  the 
case  of  a  common  submerged  body,  the  wing  span  of  the  major  control 
surfaces  is  in  the  order  of  the  body  diameter  and  the  length  of  the  body 
in  terms  of  the  wing  span  may  be  in  the  order  of  ten.  Consequently,  the 
trailing  vortices  of  a  submerged  body  are  significantly  affected  by  the 
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flow  about  the  body  and  by  the  propulsion  system  which  may  modify  the  flow 
about  the  body  so  as  to  affect  the  evolution  of  the  trailing  vortices 
and  the  occurrence  of  vortex  breakdown.  Evidently,  the  said  differences 
stem  from  the  fact  that  the  wings  on  an  aircraft  provide  the  lift  neces¬ 
sary  for  flight  whereas  the  control  surfaces  on  a  submerged  body  provide 
only  partial  lift  for  depth  and  attitude  control.  The  major  portion  of 
the  lift  is  already  provided  by  the  buoyancy  of  the  body. 

Numerous  attempts  have  been  made  to  model  the  behavior  of  laminar  and 
turbulent  trailing  vortices.  Direct  calculation  of  the  structure  of  the 
trailing  vortices  is  done  by  assuming  a  final  flow  configuration  and  re¬ 
quiring  conservation  of  some  chosen  flow  quantities.  Most  early  discus¬ 
sions  of  the  far-field  structure  were  based  on  the  vortex-wake  model  of 
Spreiter  and  Sacks  [Ref.  38].  They  have  assumed  that  the  wake  consists 
of  two  Rankine  vortices,  with  uniform  vorticity  inside  the  core  and  zero 
vorticity  outside.  The  circulation  and  vortex  separation  are  given  by  • 
Stokes'  theorem  and  conservation  of  momentum.  The  core  radius  is  deter¬ 
mined  by  requiring  kinetic  energy  to  be  conserved.  For  an  elliptically 
loaded  wing  this  model  gives  a  core  radius  of  about  0.1  times  the  trail¬ 
ing  vortex  separation.  Test  flight  data  reported  that  the  cores  are 
somewhat  smaller  than  predicted  and  that  all  of  the  vorticity  is  not  in 
the  core  [Refs.  1-2]. 

Donaldson  [Ref.  39]  noted  that  the  model  of  Betz  [Ref.  5]  which  assumes 
that  each  half  of  the  vortex  sheet  rolls  up  into  a  circular  vortex  core 
conserving  circulation,  impulse,  and  momentum  of  impulse,  could  be  used. 
This  conserves  the  moment  of  inertia  of  the  vortex  distribution  about  its 
centroid.  Betz  further  determined  the  distribution  of  circulation  in  the 
vortex  by  assuming  the  vorticity  outboard  of  some  point,  which  will  form 
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the  core  of  the  vortex  Inside  some  radius,  conserves  its  moment  of  iner¬ 
tia  about  its  own  centroid  during  roll-up.  Both  Jordan  [Ref.  40]  and 
Donaldson  [Ref.  41]  have  obtained  similar  solutions.  This  approach 
agrees  much  better  with  flight  test  data  than  the  model  of  Spreiter  and 
Sacks  [Ref.  38]. 

C.  STABILITY  OF  TRAILING  VORTICES  AND  THE  EFFECT  OF  STRATIFICATION 

If  the  trailing  vortices  could  be  represented  as  a  two-dimensional 
vortex  pair,  the  study  of  their  motion  may  be  more  straightforward. 
Unfortunately,  they  are  most  unstable  to  three-dimensional  disturbances 
so  that  the  vortices  rapidly  lose  their  two-dimensional  character.  This 
is  one  of  the  major  differences  between  the  motion  of  a  single  vortex 
and  the  vortex  pair.  A  single  vortex  is  known  to  be  very  stable  whereas 
a  vortex  pair  suffers  from  numerous  instabilities.  Nevertheless,  the 
study  of  two-dimensional  vortex  pairs  has  continued  in  hopes  of  obtaining 
a  satisfactory  solution  that  will  match  what  little  three-dimensional 
data  are  available.  The  majority  of  the  data  have  been  obtained  in  large 
scale  atmospheric  tests  using  aircraft  [Refs.  42-43]  or  from  vortex  rings 
[Refs.  44-45].  Barker  and  Crow  [Ref.  46]  and  Tomassian  [Ref.  47]  have 
created  two-dimensional  vortices  in  a  small  tank  with  large  circulations 
and  comparably  small  separations.  There  are  significant  differences 
between  the  experimental  results  cited.  These  differences  and  their 
explanations  are  further  compounded  by  the  contradictions  of  the  analyses. 

There  are  fundamental  differences  in  the  assumptions  made  and  widely 
differing  conclusions  have  been  reached  regarding  the  motion  of  the  trail¬ 
ing  vortex  pair  in  a  stratified  medium.  The  lack  of  consistent  laboratory 
data  did  not  help  to  resolve  the  differences  and  to  delineate  the  most 
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important  elements  of  the  physics  of  the  motion.  Most  models  use  an  oval 
shaped  fluid  body  carried  by  the  vortices  under  the  influence  of  their 
mutual  induction  and  the  density  difference  with  the  surroundings .  There 
is  a  significant  difference  in  the  accounting  of  fluid  momentum,  impulse, 
energy,  viscous  and  turbulent  diffusion,  and  in  the  assumptions  concern¬ 
ing  the  buoyancy  induced  vorticity  entrainment  and  detrainment  within 
the  recirculation  cell  (see  Fig.  1).  One  school  of  thought  holds  that  the 
recirculation  cell  buoys  upward  in  a  stably  stratified  atmosphere  [Refs. 
48-49],  while  others  maintain  that  the  vortices  accelerate  downward 
[Refs.  50-51].  The  resolution  of  the  dilenma  is  important,  because  of 
the  understanding  of  wake  hazard  and  of  the  energy  transferred  to  the 
surrounding  medium. 

Table  I,  adapted  from  Widnall  [Ref.  52],  describes  the  differences  in 
some  of  these  models.  The  models  of  Scorer  and  Davenport  [Ref.  50]  and 
Crow  [Ref.  51]  do  not  include  turbulence  effects.  The  models  of  Hecht 
et  al . ,  [Refs.  44  and  53]  are  computer  generated  flow  fields  assuming 
constant  recirculation  cell  areas.  Clearly,  the  differences  in  these 
models  cannot  be  resolved  until  reliable  laboratory  data  are  available. 

It  is  important  that  the  experiments  model  the  trailing  vortices  and  not 
just  a  pair  of  two-dimensional  vortices.  The  three-dimensional  nature 
of  the  flow  and  the  importance  of  the  axial  velocity  in  the  trailing 
vortices  must  be  reflected  in  the  proper  simulation  of  the  behavior  of 
vortices.  It  is,  thus,  important  that  the  differences  between  the  three- 
dimensional  experiments  and  the  two-dimensional  analytical  or  numerical 
models  be  clearly  delineated.  It  is  only  through  such  comparisons  that 
one  would  be  able  to  determine  as  to  whether  enough  physics  have  been 
incorporated  into  the  analytical  model. 


26 


Figure  1.  Downward  Migra' 
in  a  Stratifies 


1 


TABLE  I 


Comparison  of  Theoretical  Studies  of  Vortex 
Pairs  in  a  Stably  Stratified  Atmosphere 


Author 

Vortex 

Spacinq 

Migration 

Speed 

Buoyancy 

Generated 

Vorticity 

Costen  [52] 

Weak  stability 

Decreases 

Increases 

No  Effect 

Strong  stability 

Increases 

Stops 

- 

Kuhn  &  Nielsen  [58] 

Decreases 

Increases 

Partly  Ent. 

Saffman  [49] 

Increases 

Stops 

Entrained 

Scorer  &  Davenport  [50] 

Decreases 

Increases 

Detrained 

Tombach  [59] 

Strong  stability 

Increases 

Stops 

Entrained 

Weak  stability 

Decreases 

Increases 

Entrained 

Tulin  &  Shwartz  [48] 

Increases 

Stops 

Entrained  & 
Annihilated 

Crow  [51] 

Decreases 

Increases 

Created  in 
Wake 

Narain  &  Uberoi  [60] 

Decreases 

Decreases 

Entrai ned 

Hecht  et  al .  [44] 

Decreases 

Stops 

Created  at 
Boundary 

Hecht,  Bilanin, 

Hirsh  [53] 

Not 

Appreciable 

Stops 

Created  at 
Boundary 
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The  first  quantitative  analysis  of  the  three-dimensional  sinusoidal 
instability  observed  in  trailing  vortices  in  an  ideal  homogeneous  fluid 
was  given  by  Crow  [Ref.  54].  His  results  showed  good  agreement  with  the 
observations  of  the  general  features  of  aircraft  vortex  instability. 

Crow  and  Bate  [Ref.  55]  explored  the  growth  of  vortex-pair  instability 
when  driven  by  ambient  turbulence.  They  have  calculated  the  time  to 
1inkag°  as  a  function  of  the  Kolmogorov  turbulence  parameter.  The  possi¬ 
bility  of  exciting  the  vortex-pair  instability  by  changing  the  lifting 
surface  span-load  distribution  was  suggested  by  Crow  and  demonstrated  in 
flight  by  Chevalier  [Ref.  56]  by  changing  aircraft  pitch  angle  periodically. 

The  trailing  vortices  may  also  suffer  from  instabilities  inherent 
to  single  vortices,  i.e.,  from  vortex  breakdown.  Although  the  mechanism 
is  not  well  understood,  the  presence  of  axial  velocities  in  the  trailing 
vortices  and  the  proximity  of  the  body  from  which  the  vortices  are 
emanating  play  significant  roles  in  the  occurrence  of  vortex  breakdown. 

The  vortex  breakdown  has  been  studied  in  detail  for  single  vortices 
(see  e.g.,  Refs.  6-8).  It  is  not  known  what  would  be  the  effect  of  mutual 
induction  and  straining  of  vortices  on  the  breaking  of  one  or  both  of 
the  vortices.  Furthermore,  it  is  not  clear  as  to  what  happens  after  the 
breakdown.  It  was  assumed  that  the  breakdown  leads  to  the  total  dissi¬ 
pation  of  circulation  and  thus  its  occurrence  is  most  desirable.  However, 
observations  [Ref.  2]  have  shown  that  the  vortex  breakdown  does  not 
totally  dissipate  the  existing  circulation.  Thus,  the  determination  of 
how  much  circulation  remains  in  a  vortex  after  the  breakdown  constitutes 
one  of  the  major  questions.  One  would  also  like  to  determine  as  to  how 
long  the  remaining  circulation  persists.  Qualitatively,  it  is  now  clear 
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that  the  path  to  the  ultimate  destruction  of  trailing  vortices  has  not 
yet  been  charted.  The  Crow  instability,  laminar  and  turbulent  diffusion, 
entrainment  and  detrainment  of  oppositely-signed  vorticity,  vortex  break¬ 
down  or  a  cascade  of  vortex  breakdowns,  and  the  proximity  effects  (body 
or  ground)  play  significant  roles.  For  a  submerged  body,  the  deceleration 
of  the  ambient  flow  toward  the  tail  of  the  body  may  give  rise  to  vortex 
breakdown  at  a  fixed  point  relative  to  the  body.  This,  in  turn,  results 
in  a  series  of  vortex  breakdowns  as  the  body  moves  at  a  constant  speed, 
leaving  behind  a  cascade  of  breakdowns.  Since  each  breakdown  separates 
a  supercritical  swirling  flow  from  a  subcritical  flow  (downstream  of  the 
breakdown)  (see  Ref.  7)  and  since  the  subcritical  flow  cannot  sustain 
waves,  each  breakdown  rapidly  moves  upstream  to  catch  up  with  the  newest 
breakdown.  This,  in  turn,  leads  to  rapid  weakening  of  the  trailing  vortices. 
The  quantification  of  these  qualitative  observations  constitute  the  essence 
of  future  research  efforts. 

D.  SUMMARY  OF  THE  PREVIOUS  INVESTIGATIONS 

It  is  clear  from  the  foregoing  that  a  great  deal  of  additional  work 
is  needed  to  understand  the  rolling-up  of  vortices  in  both  the  strati¬ 
fied  and  unstratified  fluid  medium,  the  evolution  of  trailing  vortices, 
the  demise  mechanisms,  and  the  processes  leading  to  the  ultimate  destruc¬ 
tion  of  vortices.  Analytical  studies  must  be  based  on  reliable  experimental 
data  in  order  to  incorporate  as  much  relevant  physics  as  possible  into 
the  models.  Furthermore,  the  analyses  must  be  based  on  the  observations 
and  measurements  of  three-dimensional  trailing  vortices  rather  than  those 
of  two-dimensional  vortices  in  order  to  bring  into  focus  the  important 
role  played  by  the  three-dimensionality  of  the  flow  and  the  axial  velocity 
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in  the  vortices.  In  fact,  a  rigorous  investigation  of  the  entire  phenomenon 
might  begin  with  trailing  vortices  generated  by  relatively  simple  bodies 
(such  as  a  Delta  wing  or  rectangular  wing)  in  an  unstratified  medium  and 
then  expanded  to  include  the  effect  of  stratification.  At  present  there 
are  no  extensive  experiments  in  either  medium  and  no  detailed  calcula¬ 
tions  have  been  reported. 
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III.  EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 

It  is  clear  from  the  foregoing  that  there  is  very  little  or  no  data 
regarding  the  motion  of  trailing  vortices  in  stratified  and  unstratified 
fluids.  It  is  also  clear  that  the  analytical  models  and  numerical  calcu¬ 
lations  can  be  judged,  validated,  and  upgraded  only  through  the  use  of 
data  obtained  with  three-dimensional  trailing  vortices  under  controlled 
laboratory  conditions.  This,  in  turn,  requires  the  use  of  a  large  water 
basin  in  which  various  models  may  be  towed  at  desired  speeds. 

The  test  facility  used  in  the  present  investigation  was  developed 
with  the  following  specific  objectives  in  mind:  (i)  the  establishment  of 
the  parameters  governing  the  phenomenon  and  the  scaling  laws;  ( i i )  the 
understanding  of  the  characteristics  of  the  sinusoidal  instability  in  three- 
dimensional  trailing  vortices  with  axial  velocity;  (iii)  the  determination 
of  the  role  played  by  the  vortex  breakdown  or  cascade  of  breakdowns  and  of 
the  circulation  remaining  in  the  trailing  vortices  following  breakdown; 

(iv)  the  establishment  of  any  interaction  between  the  sinusoidal  instability 
and  the  vortex  breakdown;  (v)  the  investigation  of  the  effect  of  ambient 
stratification,  ambient  turbulence,  and  model-induced  turbulence  on  the 
evolution  and  demise  of  trailing  vortices;  (vi)  the  determination  of  the 
characteristics  of  internal  waves  generated  by  the  body  and/or  trailing 
vortices  prior  to,  during,  and  after  the  onset  of  sinusoidal  instability 
and/or  vortex  breakdown;  (vii)  to  examine  the  behavior  of  trailing  vortices 
generated  by  a  lifting  surface  undergoing  time-dependent  oscillations  (with 
or  without  the  body-proximity  effects);  (vlli)  to  study  the  evolution  of 
Horseshoe  vortices  and  vortex  rings;  (lx)  to  examine  the  effect  of  the 
alterations  in  the  axial  velocity  of  the  vortices  at  the  time  of  their 
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generation;  and  (x)  to  investigate  the  interaction  of  the  trailing  vortices 
with  the  vortices  generated  by  the  body  and/or  control  surfaces. 

The  experimental  equipment  consists  of  the  water  basin,  the  strati¬ 
fication  system,  and  the  test  models. 

A.  WATER  BASIN 

The  basin  is  36  feet  long,  3  feet  wide  and  5  feet  deep  and  is  made  of 
0.5-inch  thick  aluminum  plates,  welded  appropriately  (see  Fig.  2).  The 
back  wall  of  the  tank  is  reinforced  to  prevent  bulging  under  hydrostatic 
pressure.  The  front  wall  Is  fitted  with  square  windows  of  1.0-inch  thick 
plexiglass  for  flow  visualization. 

Vertical  and  horizontal  scales  of  0.5-inch  increments  are  marked  on 
the  central  window.  The  scales  are  used  during  the  experiments  to  evaluate 
the  position  of  the  trailing  vortices  (see  Fig.  3). 

The  aluminum  interior  of  the  tank  is  painted  with  marine  epo)<y  to 
prevent  corrosion.  In  addition,  the  inside  face  of  the  back  wall  is  painted 
with  flat  white  paint  to  provide  better  contrast  and  make  the  flow  visuali¬ 
zation  easier  with  the  dye  used. 

The  filling  pipes  lie  on  the  bottom  of  the  tank.  They  are  plastic 
pipes  2.0-inch  diameter,  perforated  along  their  entire  length  with  holes 
of  1/16-inch  diameter  (see  Fig.  4). 

Two  drains  are  provided  at  the  bottom  at  each  end  of  the  basin.  Two 
parallel  rails  of  1.0-inch  diameter  are  mounted  along  the  bottom  of  the 
tank.  A  carriage  rides  smoothly  on  these  rails  and  provides  the  test 
body  with  a  constant  velocity  through  the  use  of  an  endless  cable  and  a 
3/4  Hp  DC  motor  (see  Fig.  5).  The  velocity  of  the  model  is  measured  and 
continuously  monitored  through  the  use  of  a  magnetic  velocity  sensor. 
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Figure  2.  Water  Basin 


iyure  3.  Central  Window  and  the  Horizontal  and  Vertical  Scales. 


Figure  4.  Top  View  of  Strati  fixation  System. 


mounted  above  a  precision  gear.  The  sensor  counts  the  number  of  teeth 
and  yields  and  instantaneous  velocity  proportional  to  the  velocity  of 
the  model . 

The  two  rails,  the  carriage  and  the  filling  pipes  are  located  on  or 
near  the  bottom  of  the  basin  and  under  a  turbulence  management  system. 

This  system  consists  of  a  one-inch-thick  polyurethane  foam,  sandwiched 
between  two  perforated  aluminum  plates  (see  Fig.  6). 

B.  STRATIFICATION  SYSTEM 

The  stratification  system  consists  of  four  tanks  (60-gallon  capacity 
each)  and  the  necessary  tubes  and  valves. 

First,  brine  water  of  certain  salt  concentration  is  prepared  in  one 
of  these  tanks.  Then  precalculated  amounts  of  brine  water  are  added  to 
the  other  three  tanks,  depending  on  the  desired  density  gradient  and  the 
number  of  the  layers.  The  calculation  of  the  amount  of  salt  needed,  for 
a  certain  density  gradient  or  for  a  certain  stratification  parameter,  is 
discussed  in  Appendix  A. 

The  resulting  solution  is  drained  into  the  basin  at  a  rate  of  2.1  CFM. 
Necessary  precautions  have  been  taken  to  prevent  air  entrainment  and 
subsequent  bubble  generation.  Furthermore,  two  air-bubble  traps  have 
been,  installed,  one  just  before  the  ma’’n  valve  of  the  filling  system 
and  one  at  the  middle  point  of  the  delivery  pipes. 

This  system  can  fill  the  water  basin  and  establish  a  desired  stratifi¬ 
cation  in  about  three  hours.  There  is  no  difficulty  regarding  the  desired 
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density  gradient,  even  for  the  lower  gradients  (e.g.,  do/dy  =  5x10 
slug/ ft  /inch).  The  density  discontinuity  between  the  0.9- inch  layers 
disappeared  by  diffusion  after  a  short  time  for  the  stratifications  of 
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interest.  Smooth  density  gradients  were  obtained  as  verified  by  a 
salinity  meter.  The  parts  of  the  stratification  system  are  shown  in 
Figs.  4  and  7. 

C.  MODELS 

Three  models  were  used  in  the  experiments.  The  first  two  models  are 
Delta  wings,  made  of  1/8-inch  thick  aluminum  with  hollow  interiors. 

Each  model  cavity  is  filled  with  dye  prior  to  the  filling  of  the  basin 
with  water.  The  two  holes  at  the  rear  edge  of  each  model  and  the  third 
hole  in  the  middle  of  the  upper  surface  are  plugged  prior  to  an  experi¬ 
mental  run  to  prevent  dye  leakage  (see  Fig.  8).  The  dye  consisted  of 
food  coloring  and  proper  amounts  of  alcohol.  The  density  of  the  mixture 
was  such  that  the  dye  was  only  extremely  slightly  buoyant  relative  to 
its  surroundings. 

The  models  are  mounted  on  their  bases  by  means  of  a  streamlined  thin 
aluminar  bar  (its  cross  section  was  a  NACA  0006  foil)  and  set  at  the  de¬ 
sired  angle  of  attack.  The  models  are  placed  at  a  level  9  inches  above 
the  foam.  The  first  Delta  wing  had  a  vertex  angle  of  39.5055  degrees  and 
a  base  length  of  7.9  inches.  The  second  model  had  a  vertex  angle  of  50 
degrees  and  a  base  length  of  5.071  inches. 

The  third  model  is  a  submerged  planing  surface  made  of  plexiglass  and 
attached  to  a  streamlined  body.  The  body  is  cut  in  half  longitudinally 
with  the  bottom  half  removed  to  allow  testing  at  greater  depths  with 
negligible  hydrodynamic  ground  effects.  The  planing  surface  has  a  base 
length  of  6.96  inches  with  an  aspect  ratio  of  2.20  (Aspect  ratio  =  [base 
length]  /Area).  The  model  is  mounted  to  the  carriage  such  that  the  plane 
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Figure  8.  Delta-Wing  Models. 


bottom  of  the  body  is  flush  with  the  top  of  the  aluminum  plate  above  the 
foam. 

A  hollowed  core  within  the  supporting  column  of  the  planing  surfaces 
is  filled  with  dye.  The  dye  is  dispensed  from  the  tips  of  the  planing 
surfaces  with  small  holes  connected  to  the  dye  reservoir  (i.e.,  the 
hollow  core) . 

As  noted  earlier  all  models  are  pulled  by  means  of  a  DC  motor,  pulley, 
and  cable  system  at  the  desired  speed  (ranging  from  1.0  to  5.0  ft/sec). 

The  entire  arrangement  is  shown  in  Figs.  9  and  10. 

0.  TEST  PROCEDURES 

The  model  is  initially  filled  with  dye  and  all  the  fill  and  drain 
holes  on  the  model  are  closed  with  wooden  pegs.  Then  the  basin  is  filled 
gradually  either  with  fresh  water  (for  unstratified  flow  experiments)  or 
with  salt  water  (for  stratified  flow  experiments)  to  the  desired  level. 
After  sufficient  period  of  waiting  for  the  elimination  of  any  internal 
currents  in  the  basin  the  model  was  set  in  motion  at  a  desired  speed. 

The  motion  of  the  trailing  vortices  are  recorded  on  high-speed  Polaroid 
film  at  the  test  section  (one  of  the  plexiglass  panels  near  the  middle 
of  the  basin).  Each  picture  included  two  clocks  accurate  to  0.01  second, 
the  vertical  and  horizontal  scales  on  the  plexiglass  window  and,  of 
course,  the  side  view  of  the  trailing  vortices  as  they  rose  from  the  model 
after  formation  (see  Fig.  11).  The  time  interval  between  successive  pic¬ 
tures  is  determined  from  the  two  clocks.  The  vertical  rise  of  the  vortices 
is  determined  from  the  vertical  scale.  Attention  has  been  paid  to  the 
fact  that  the  vortices  are  farther  away  from  the  scale  on  the  plexiglass 
window  and  that  the  scale  placed  vertically  in  the  middle  of  the  test 
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Figure  9.  Motion  Generation  Mechanism. 


section  does  not  exactly  correspond  to  the  scale  marked  on  the  plexi¬ 
glass  window  due  to  refraction  and  parallax.  The  necessary  correction 
was  made  by  photographing  a  scale  placed  in  water  in  the  middle  of  the 
test  section  together  with  the  scale  marked  on  the  plexiglass  window. 

This  resulted  in  a  simple  conversion  table  which  enabled  the  determina¬ 
tion  of  the  actual  position  of  the  vortex  from  the  scale  reading  on  the 
photograph. 

The  results  are  non-dimensional ized  and  plotted  in  various  forms  and 
compared  with  those  obtained  in  the  previous  runs.  Each  experiment  was 
repeated  at  least  twice  to  ascertain  that  the  experiments  were  repeatable. 

The  experiments  with  the  formation  of  Horseshoe  vortices  and  vortex 
rings  were  conducted  in  the  manner  described  above  with  the  exception  that 
the  model  was  initially  placed  closer  to  the  test  section,  rather  than 
at  one  end  of  the  basin,  and  that  both  the  horizontal  and  vertical  scales 
were  used  to  determine  the  velocity  components  of  the  translation  of  the 
vortices. 

All  trailing  vortices  and  vortex  rings  were  recorded  on  film  until  the 
time  they  have  completely  dissipated  either  due  to  aging  (diffusion  of 
vorticity  due  to  viscosity,  turbulence,  entrainment,  and  detrainment) ,  or 
due  to  Crow  instability  (sinusoidal  instability  and  linking  leading  to 
the  formation  of  vortex  rings),  and/or  due  to  vortex  breakdown  (core 
bursting).  It  was  thus  possible  to  determine  the  life  span  of  vortex 
cores  from  their  formation  to  their  final  demise. 
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IV.  RESULTS 


A.  DIMENSIONAL  ANALYSIS 

The  experimental  procedure  is  based  on  a  careful  dimensional  analy¬ 
sis  of  the  parameters  entering  into  the  phenomenon. 

The  dependent  parameter  of  major  importance  is  the  displacement  H 
of  the  vortex  in  the  vertical  direction  and  it  may  be  expressed  as  a 
function  of  the  following  independent  parameters  for  a  given  body  shape: 


H  =  f(t,U,D,p0,|£-,v,B,e,a,g) 


(1) 


where  t  denotes  the  time;  U,  the  forward  velocity  of  the  model;  D,  the 
initial  depth  of  submergence  of  the  vortex  pair;  pQ,  the  fluid  density; 
3p/3y,  the  density  gradient;  \j,  the  kinematic  viscosity  of  water;  B,  the 
base  width  of  the  model  (see  Fig.  8);  e,  the  vertex  angle;  a,  the  angle  of 
attack;  and  g,  the  gravitational  acceleration.  The  other  parameters  such 
as  the  vortex  strength  r*  the  core  radius  of  the  vortices  rQ,  the  vortex 
separation  bQ,  and  the  velocity  of  vertical  migration  of  vortices  VQ, 
are  determined  by  the  same  parameters  governing  the  variation  of  H. 

A  dimensional  analysis  of  Eq.  (1)  yields 
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and  is  known  as  the  Brunt-Vaisala  frequency.  The  other  dimensionless 
parameters  appearing  in  Eq.  (2)  are  self-explanatory. 


In  the  present  investigation  B,  e,  a,  and  D  were  determined  by  the 
model  being  investigated.  The  velocity  of  the  model  was  varied  from 
0.32  ft/s  to  4.04  ft/s  in  experiments  with  fresh  water  (N  =  0).  The 
experiments  with  stratified  water  were  carried  out  by  varying  N 
(N  =  0.0974  s"1  to  0.2101  s’^)  at  various  velocities. 

The  dimensionless  parameter  NB/U  is  often  replaced  by  another  parameter 
NbQ/Vo>  known  as  the  stratification  parameter.  It  should  be  noted  that 
both  bQ  and  VQ  are  dependent  parameters.  bQ  is  determined  by  the  shape  of 
the  body  and  its  angle  of  attack.  VQ  is  determined  by 
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Thus,  Eq.  (2)  may  be  recast  as 
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Of  the  parameters  appearing  in  Eq.  (4)  only  t,  D,  N,  g,  and  v  are  inde¬ 
pendent  and  may  be  varied  irrespective  of  the  shape  of  the  body.  However, 
the  remaining  parameters  with  the  exception  of  a  depend  on  the  shape, 
aspect  ratio,  and  the  velocity  of  the  body.  The  variation  of  e  and  a 
are  imbedded  in  the  variation  of  VQ  and  b  .  Experiments  were  performed 
with  a  given  body  by  varying  the  angle  of  attack  and  the  speed  of  the 
body  for  a  given  Brunt-Vaisala  frequency,  N.  Then  VQ  and  H  were  deter¬ 
mined  from  the  pictures.  The  parameters  e  and  bQ  were  obtained  from  the 


existing  experimental  data  [Ref.  61].  Subsequently,  the  data  were  plotted 

in  terms  of  H/bQ  versus  VQt/bo  for  various  values  of  D/bQ,  and  Nbo/VQ. 

2 

The  measurements  have  shown  that  the  dimensionless  parameters  Vo/gbo  and 
VQb0/v  are  not  important  within  the  range  of  the  parameters  encountered 
in  the  experiments.  With  the  foregoing  arguments,  Eq.  (4)  may  be  written 
as 


U  V  t  n  Nbn 

H  fi  0  _D _ 0\ 

b  '  h  ’b  ’  V  ; 

o  ooo 


(5) 


This  equation  formed  the  basis  for  the  analysis  of  the  entire  experimental 
data. 


B.  DELTA-WING  GENERATED  VORTICES  IN  UNSTRATI FIED  WATER 

Sample  data  obtained  with  the  first  Delta-wing  (referred  to  hereafter 
as  Delta  1)  in  unstratified  water  (NbQ/VQ  *  SP  =  0)  are  shown  in  Fig.  12 
for  two  values  of  VQ  and  for  D/bQ  *  5.19.  Two  symbols,  identifying  two 
independent  runs,  show  that  there  is  very  little  scatter  in  the  data.  Had 
the  vortices  maintained  their  strength  and  spacing,  as  they  would  in  an 
ideal  fluid,  H/bQ  would  have  been  equal  to  VQt/bo .  This  straight  line 
relationship  is  also  shown  in  Fig.  12.  Clearly,  all  the  data  points  are 
below  the  ideal  line  and  show  that  the  vortices  do  not  retain  their  original 
characteristics.  One  may  argue  that  this  may  be  as  a  consequence  of 
several  factors  such  as  the  decrease  of  circulation,  the  increase  of  the 
vortex  spacing  b,  and  the  proximity  of  the  free  surface.  The  variation  of 
b  with  time  may  be  excluded  from  this  discussion  because  the  observations 
and  measurements  have  shown  that  it  remains  practically  constant  for  large 
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Figure  12.  Vertical  Vortex  Rise  for  Delta  1,  SP  =  0,  D/b  =  5.19. 


enough  values  of  D/bQ.  As  to  the  free  surface  effects,  it  will  be  shown 
that  they  do  not  play  any  role  in  the  rise  of  the  vortices  for  0/bQ 
larger  than  about  4.3.  Thus  one  may  conclude  that  for  large  D/bQ  values, 
the  difference  between  the  actual  and  the  ideal  values  of  H/bQ  is  due  to 
the  decrease  of  the  circulation  of  vortices  with  time.  This  decrease,  as 
noted  in  connection  with  the  discussion  of  previous  studies,  is  a  conse¬ 
quence  of  eddy  diffusion  of  vorticity,  Crow  instability,  and  vortex  break¬ 
down.  This  fact  has  not  been  recognized  in  the  past  because  of  lack  of 
reliable  data  and  it  was  often  assumed  that  the  vortex  strength  remained 
constant.  This  in  turn  gave  rise  to  conflicting  conclusions  cited  in 
Table  I. 

Figure  13  shows  a  composite  plot  of  the  data  obtained  with  Delta  1 
for  various  values  of  VQ  at  a  fixed  value  of  D/bQ  *  4.31.  The  primary 
purpose  of  this  figure  is  to  substantiate  the  assertion  that  the  dimension 
less  parameters  VQ/gbo  (a  representative  Froude  number)  and  VQbo/v  (a 
representative  Reynolds  Number)  do  not  play  experimentally  detectable 
roles  in  the  variations  of  H/bQ  with  VQt/bo>  at  least  for  unstratified 
flow  experiments.  The  two-fold  variation  of  VQ  for  constant  values  of  bQ, 
g,  and  v  yields  a  four-fold  variation  in  the  Froude  number  and  a  two-fold 
variation  in  the  Reynolds  number.  Clearly,  Fig.  13  does  not  exhibit  such 
a  dependence  on  the  said  two  parameters.  In  fact,  it  is  because  of  this 
finding  that  the  governing  parameters  were  finally  reduced  to  that  given 
by  Eq.  (5). 

The  next  question  to  be  addressed  is  the  dependence  of  H/bQ  on  D/bQ. 

In  other  words,  at  what  0/bQ  does  the  free  surface  affect  the  rise  of  the 
vortices  and  in  what  form  does  this  effect  exhibit  itself? 
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Figure  13.  Vertical  Rise  for  Delta  1,  SP  =  0,  D/b  =  4.31 


Figure  14  shows  the  data  obtained  with  0/bQ  values  ranging  from  2.15 
to  5.19,  again  in  unstratified  water.  Evidently,  the  variation  of  H/bQ 
versus  VQt/b0  for  D/bQ  =  4.31  and  D/bQ  =  5.19  is  such  that  the  differ¬ 
ences  are  negligible  and  certainly  well  within  the  experimental  errors. 

For  other  values  of  D/bQ,  the  vortices  reach  their  limiting  height  very 
quickly  and  then  begin  to  bounce  back  from  the  free  surface.  The  smaller 
the  0/bQ,  the  smaller  is  the  limiting  height.  These  experimental  facts 
may  be  explained  partly  in  terms  of  the  behavior  of  inviscid  vortices  in 
the  vicinity  of  a  plane  wall. 

The  free  surface  acts  as  a  reflection  plane  and  the  motion  of  the 
vortex  pair  is  determined  not  by  their  mutual  induction  alone  but  also  by 
the  image  vortex  pair  above  the  free  surface.  In  other  words,  the  tra¬ 
jectory  of  the  ideal  vortices  is  determined  by  the  mutual  induction  of 
four  vortices  (see  Fig.  15).  Lamb  [Ref.  62]  has  shown  that  the  trajectories 
of  line  vortices  in  ideal  unstratified  fluid  are  hyperbolas  given  by 


(2x/bo)‘ 


(2y/bo)‘ 


=  1 


.  f bO\2 

+  {7D} 


(6) 


where  bQ  is  the  initial  spacing  of  the  two  real  vortices,  and  D  is  their 
initial  depth.  Equation  (6)  shows  that  the  vortices  move  sideways  as 
they  approach  the  free  surface.  For  large  values  of  x  and  D,  y  approaches 
half  the  initial  separation,  i.e.,  the  vortex  pair  does  not  come  closer 
to  the  free  surface  than  a  distance  of  bQ/2 . 

A  similar  analysis  cannot  yet  be  carried  out  in  a  real  fluid  including 
the  effects  of  stratification  and  turbulence.  The  strength  of  the  vortices 
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Figure  14.  D/b  Effects  on  Vertical  Vortex  Rise  for  Delta 


Figure  15.  Mutual  Induction  of  Four  Vortices  and  the 
Free  Surface  Effects. 
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and  their  rise  velocity  continue  to  decrease  (see  the  slope  of  the  data 
for  0/bQ  *  2.15  in  Fig.  14). 

As  the  vortices  move  too  far  sideways,  they  begin  to  fall  under  the 
influence  of  additional  image  vortices  due  to  the  presence  of  the  side 
walls.  Even  though  it  is  not  a  matter  of  primary  importance,  one  can 
show  that  the  combined  effect  of  the  image  vortices  due  to  the  side  walls 
and  the  free  surface  is  to  push  the  vortices  downward.  It  is  not  yet 
possible  to  verify  whether  there  exists  a  physical  mechanism  which  will 
make  the  vortices  bounce  from  the  free  surface  had  there  been  no  side 
wall  effects.  It  is  possible  that  the  decrease  in  the  mutual  influence 
of  the  trailing  vortices  may  give  rise  to  alterations  in  the  velocity 
fields  in  the  presence  of  a  movable  free-surface  boundary  and  this  in 
turn  may  result  in  a  bouncing  back  of  the  vortices.  In  any  case,  only 
the  experiments  in  a  large  basin  can  resolve  the  question. 

An  additional  fact  which  can  be  extracted  from  Fig.  14  is  the  lack 
of  importance  of  the  Froude  and  Reynolds  numbers  on  the  rise  of  vortices 
even  when  the  free  surface  effects  predominate.  The  data  for  all  D/bQ 
values  were  obtained  at  different  model  speeds  even  though  they  have  been 
shown  with  a  single  symbol  in  Fig.  14.  The  plotted  as  well  as  the  tabu¬ 
lated  data  (see  Appendix  C)  show  that  (see  for  example  the  data  for 
D/bQ  *  2.15  in  Fig.  14)  neither  the  Froude  number  nor  the  Reynolds  number 
plays  any  role  on  the  rise  of  vortices  in  unstratified  flow  with  or 
without  free  surface  effects,  within  the  range  of  the  parameters  encoun¬ 
tered  In  this  investigation. 

The  results  presented  above  were  obtained  with  Delta  1.  In  order  to 
ascertain  further  the  validity  of  the  conclusions  reached  regarding  the 
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role  played  by  the  governing  parameters  (see  Eq.  (5))  experiments  were 
carried  out  with  the  second  Delta-wing  model  (Delta  2).  Figure  16  shows 
H/bQ  versus  VQt/b0  for  D/bQ  =  7.73  in  unstratified  water  for  various 
values  of  VQ  (Tabulated  data  are  presented  in  Appendix  D).  Clearly,  the 
results  do  not  depend  on  VQ  and  hence  on  the  model  speed.  The  comparison 
of  the  data  obtained  with  the  two  models  is  presented  in  Fig.  17.  Aside 
from  some  minor  scatter  at  large  values  of  V0t/bQ,  the  results  of  the 
two  models  compare  quite  well.  Thus,  it  may  be  concluded  on  the  basis 
of  the  extensive  data  collected  from  the  two  models  that  in  unstratified 
water  the  rise  of  the  trailing  vortices  is  governed  by  the  dimensionless 
parameters  H/bQ,  V0t/bQ,  and  D/bQ.  In  all  cases  the  vortices  gradually 
lose  their  strength  and  thus  they  move  upwards  with  ever  decreasing 
velocities  until  they  are  dissipated  by  various  demise  mechanisms  dis¬ 
cussed  earlier. 

C.  DELTA-WING  GENERATED  VORTICES  IN  STRATIFIED  WATER 

Experiments  in  stratified  water  were  carried  out  with  the  same  Delta¬ 
wing  models  for  various  stratification  parameters  (SP  =  NbQ/VQ).  First, 
the  results  obtained  with  Delta  1  for  five  SP  values  will  be  presented. 
Then  the  results  obtained  with  Delta  2  for  three  SP  values  will  be 
shown.  Then  the  results  obtained  with  the  two  models  will  be  compared. 

Figures  18  through  22  show  H/bQ  versus  VQt/bQ  for  SP  values  of  0.375, 
0.50,  0.625,  0.75,  and  1.0  for  a  fixed  value  of  D/bQ  =  4.31.  It  is 
important  to  note  that  this  particular  value  of  D/bQ  was  chosen  so  as  to 
exclude  the  effects  of  free  surface  and  thus  to  concentrate  only  on  the 
effects  of  linear  stratification.  A  complete  plot  of  the  aforementioned 
data  is  presented  in  Fig.  23.  Clearly,  the  larger  the  stratification 
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Figure  16.  Vertical  Vortex  Rise  for  Delta  2,  SP  =  0,  D/b  =  7.73. 
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Figure  17.  Comparison  of  the  Vortex  Rise  for  Delta  1  and  Delta  2  for  SP 
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Figure  18.  Vertical  Vortex  Rise  for  Delta  1,  SP  =  0.375,  D/b  =  4.31 
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Figure  19.  Vertical  Vortex  Rise  for  Delta  1 ,  SP  =  0.50,  D/b  =  4.31. 


Figure  20.  Vertical  Vortex  Rise  for  Delta  1 ,  SP  =  0„625,  D/b  =  4.31 . 
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Figure  21.  Vertical  Vortex  Rise  for  Delta  1 ,  SP  =  0.75,  D/b  =  4.31. 
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Figure  22.  Vertical  Vortex  Rise  for  Delta  1,  SP  =1.0,  D/b  =  4.31. 
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Figure  23„  Stratification  Effects  on  Vortex  Rise  for  Delta 


parameter,  the  smaller  the  rise  of  vortices.  The  last  data  points  (of 
each  particular  run)  correspond  to  the  state  where  the  vortex  core  was  no 
longer  discernible  and  the  surrounding  dye  was  randomly  diffused  with  no 
sign  of  circulation.  As  will  be  discussed  in  detail  later,  this  particu¬ 
lar  state  was  reached  after  the  vortices  were  subjected  to  one  or  more 
of  the  demise  mechanisms.  The  fundamental  conclusion  to  be  deduced  from 
Fig.  23  is  that  the  effect  of  the  stable  linear  stratification  is  to 
inhibit  the  rise  of  the  trailing  vortices.  This  may  be  due  to  a  number 
of  reasons  such  as:  (i)  the  increase  of  the  downward  buoyant  force  on 
the  recirculation  cell,  (ii)  the  generation  of  countersign  vorticity  out¬ 
side  the  recirculation  cell  (to  be  discussed  later  in  Section  IV. F), 

(iii)  the  enhanced  effect  of  the  entrainment  and/or  detrainment  of  the 
recirculation  cell,  (iv)  the  possible  enhancement  of  the  vortex  breakdown 
and  alterations  in  the  circulation  distribution  after  vortex  breakdown 
due  to  stratification,  (v)  the  change  of  the  turbulence  structure  and 
hence  the  eddy  diffusion  by  stratification,  and  (vi)  the  enhancement  of 
the  Crow  instability,  or  other  as  yet  unknown  demise  mechanisms  associated 
with  stratification. 

It  is  evident  from  Fig.  23  that  the  vortices  in  a  stratified  medium 
rise  to  smaller  heights.  Thus,  the  free  surface  effects  become  important 
at  D/bQ  values  smaller  than  that  corresponding  to  the  unstratified  case. 
Clearly,  the  particular  value  of  D/bQ  at  which  the  free-surface  effects 
become  important  in  stratified  medium  depends  on  the  stratification 
parameter. 

The  results  obtained  with  Delta  2  are  shown  in  Figs.  23  through  26  for 
three  SP  values  (SP  =  0.375,  0.50,  and  0.75)  and  for  0/bQ  =  6.50. 
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Figure  24.  Vertical  Vortex  Rise  for  Delta  2,  SP  =  0.375,  D/b  =  6.50. 
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Figure  25.  Vertical  Vortex  Rise  for  Delta  2,  SP  =  0.50,  D/bQ  =  6.50. 


Figure  26.  Vertical  Vortex  Rise  for  Delta  2,  SP  =  0.75,  D/b  =  6.50. 


Figure  25  shows  the  repeatability  of  the  data  and  all  three  figures  show 
the  general  consistency  of  the  experiments. 

Figure  27  is  a  composite  plot  of  the  aforementioned  figures.  A 
perusal  of  this  figure  reveals  the  same  conclusions  arrived  at  earlier 
in  connection  with  the  discussion  of  Delta  1  (Fig.  23). 

Since  the  purpose  of  the  use  of  two  Delta-wing  models  was  to  confirm 
the  independence  of  the  data  from  the  size  of  the  model.  Figs.  28  through 
30  were  prepared  for  comparison  of  the  data  obtained  from  Delta  1  and 
Delta  2.  In  general,  the  data  compare  extremely  well  for  a  given  strati¬ 
fication  parameter.  The  minor  differences  in  H/bQ  for  large  values  of 
V0t/b0  (see  Figs.  29  and  30)  are  thought  to  be  due  to  the  somewhat  random 
influence  of  the  demise  mechanisms  at  the  later  stages  of  the  life  of 
the  trailing  vortices.  Finally,  it  should  be  noted  that  D/bQ  does  not 
affect  the  variation  of  H/bQ  with  V0t/bQ  for  D/bQ  larger  than  at  least 
4.3,  as  discussed  earlier. 

D.  DEMISE  OF  DELTA-WING  GENERATED  TRAILING  VORTICES 

The  discussion  of  the  demise  of  the  trailing  vortices  will  be  carried 
out  first  with  regard  to  the  initial  stages  of  motion  and  then  to  the 
fully  developed  trailing  vortices. 

1 .  Horseshoe  Vortices  and  Vortex  Rings 

As  the  lifting  surface  is  set  in  motion  impulsively  from  rest  the 
vorticity  shed  from  the  surface  forms  a  so-called  Horseshoe  vortex. 
Shortly  after  the  formation  of  this  vortex,  linking  occurs  and  the 
Horseshoe  vortex  transforms  into  a  vortex  ring  (see  Fig.  31,  where  the 
Horseshoe  vortex  is  at  the  point  of  linking).  The  observations  and 
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Figure  27.  Stratification  Effects  on  Vortex  Rise  for  Delta 


Figure  28.  Comparison  of  the  Vortex  Rise  for  Delta  1  and  Delta  2  for  SP  =  0.375. 
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Figure  29.  Comparison  of  the  Vortex  Rise  for  Delta  1  and  Delta  2  for  SP  =  0.50 
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Figure  30.  Comparison  of  the  Vortex  Rise  for  Delta  1  and  Delta  2  for  SP  =  0.75. 


Figure  31.  Horseshoe  Vortex  and  the  Formation  of  a  Vortex  Ring. 


measurements  have  shown  that  the  vortex  ring  inclines  itself  at  an  angle 
of  about  45°  and  undergoes  complex  deformations  partly  under  the  influ¬ 
ence  of  self-induced  velocities  and  partly  due  to  the  velocities  created 
by  the  remainder  of  the  trailing  vortex.  The  complexity  of  the  fluid 
motion  in  the  vortex  ring  is  such  that  analyses  based  on  circular  vortex 
rings  of  constant  core  and  pitch  diameter  cannot  be  used  to  predict  the 
transport  velocity  of  the  vortex  ring.  In  fact,  attempts  made  to  do  so 
have  resulted  in  an  arbitrary  selection  of  the  core  radius  with  widely 
differing  velocities.  The  non-uniformity  of  the  velocities  induced  on  the 
ring  resulted  in  vortex  stretching  in  addition  to  gross  deformations. 

Following  the  formation  of  the  first  ring,  the  end  of  the  trailing 
vortices  necessarily  links  and  forms  a  new  Horseshoe  vortex  because  of 
the  fact  that  the  vortices  cannot  end  in  a  fluid.  This,  in 
turn,  results  in  a  new  vortex  ring  with  similar  consequences .  Repeated 
experiments  have  shown  that  as  many  as  three  or  four  such  rings  are  formed. 
Several  facts  must  be  brought  out  in  connection  with  their  formation, 
motion,  and  subsequent  demise.  Each  ring  has  a  character  of  its  own  be¬ 
cause  of  a  number  of  reasons.  The  first  ring  has  the  largest  circulation 
because  an  accelerating  body  gives  rise  to  the  largest  accumulation  of 
vorticity.  Second,  the  first  ring  is  formed  from  that  portion  of  the 
trailing  vortex  which  is  relatively  young  or  not  too  much  affected  by 
aging.  Lastly,  the  first  ring  is  under  the  influence  of  its  own  self- 
induced  velocity  plus  that  of  the  Horseshoe  vortex  at  the  end  of  the 
trailing  vortex. 

The  second  ring  is  formed  from  that  part  of  the  vortex  trail  which 
is  relatively  more  aged.  In  other  words,  the  aging  and  deformation  of 
a  vortex  ring  formed  from  a  relatively  young  vortex  trail  (higher 
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concentration  of  vorticity  at  the  core)  are  not  the  same  as  that  of  a 
ring  formed  from  an  already  aged  core.  Furthermore,  the  second  ring  is 
influenced  not  only  by  the  velocities  induced  by  itself,  but  also  by  the 
velocities  induced  by  the  first  ring  and  the  remainder  of  the  trailing 
vortex.  Extending  this  discussion  to  three  or  four  rings  one  realizes 
that  each  ring  is  different  at  birth  and  is  subjected  to  different  velocity 
fields  subsequently.  Also,  one  may  surmise  that  the  lack  of  continuation 
of  the  formation  of  a  larger  number  of  vortex  rings  may  be  related  to  the 
demise  of  the  trailing  end  of  the  vortices  by  simple  eddy  diffusion  of 
vorticity. 

From  the  foregoing,  one  may  postulate  that  the  rings  will  form 
at  such  a  frequency  that  the  total  length  of  the  trailing  vortex  between 
the  body  and  the  last  ring  will  remain  either  constant,  or  will  get  larger 
or  smaller.  Experiments  have  shown  that  not  only  is  the  number  of  rings 
finite,  but  also  the  rate  of  horizontal  progression.  The  distance  as 
well  as  the  time  interval  between  the  successive  linkings  have  been  deter¬ 
mine  and  it  was  found  that  (Ax)n/(At)n  is  nearly  constant  (Axn  is  the 
horizontal  distance  between  the  n-th  and  (n-l)-th  link  and  Atn  is  the 
corresponding  time  difference).  Comparing  the  linking  velocity  given  by 
(Ax)n/(At)n  with  the  velocity  of  the  model,  it  has  been  found  that 

wrjm ~  '  16-20  (7> 

Thus,  one  may  conclude  that  the  length  of  the  vortex  trail  increases 
rapidly  and  its  demise  is  governed  by  other  mechanisms  unrelated  to  the 
formation  of  the  Horseshoe  vortices  and  vortex  rings.  However,  the 
formation  of  the  vortex  rings  is  extremely  important  since  the  formation 
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of  the  rings  is  closely  related  to  the  change  of  circulation  generated  by 
the  lifting  body.  Such  a  change  may  result  either  from  the  change  of  the 
angle  of  attack  of  the  lifting  surface  and/or  from  the  change  of  velocity 
of  the  body.  One  may  also  include  in  the  change  of  angle  of  attack  the 
alteration  of  the  positions  of  the  flaps  of  an  aircraft.  Consequently, 
an  aircraft  or  a  submerged  body  moving  with  rapidly  changing  velocities 
or  at  a  constant  velocity  with  periodically  oscillating  lifting  surfaces 
can  give  rise  to  vortex  rings.  Finally,  it  is  worth  noting  that  the 
changes  in  velocity  and/or  the  angle  of  attack  of  lifting  surfaces  may 
precipitate  and  enhance  other  demise  mechanisms. 

2.  Demise  of  the  Vortex  Pair 

The  maximum  height  to  which  the  vortex  pair  rises  has  been  obtained 
from  the  plotted  as  well  as  the  tabulated  data  and  is  shown  in  Fig.  32  as 
a  function  of  the  stratification  parameter  Nbo/VQ.  Only  those  points  for 
which  D/bQ  were  sufficiently  large  so  as  to  preclude  any  free  surface 
effects  are  ~hown  in  Fig.  32.  Also  shown  in  this  figure  are  the  couple 
of  data  points  obtained  by  Tomassian  [Ref.  47]. 

Figure  32  points  out  several  important  facts.  First,  the  vortices 
have  a  sel f-1 imiting  effect  on  their  rise,  even  in  an  unstratified  medium. 
Second,  the  stratification  imposes  further  limitation  on  the  maximum 
height  of  rise.  In  all  cases,  various  demise  mechanisms  bring  about  the 
said  limitation.  These  will  now  be  discussed  in  some  detail,  bearing  in 
mind  the  fact  that  it  is  in  some  cases  extremely  difficult  to  quantify 
some  of  these  mechanisms. 

The  most  obvious  of  all  the  mechanisms  is  the  so-called  aging  of 
the  vortices  or  the  diffusion  of  vorticity  by  viscous  and  turbulent  actions. 


. . . 


79 


It  is  a  well-known  fact  that  a  vortex  with  a  given  initial  vorticity 
distribution  diffuses  its  vorticity  gradually.  If  the  vortices  are 
moving  in  a  stratified  medium  the  density  gradient  gives  rise  to  opposite- 
signed  vorticity  outside  the  recirculation  cell.  This  vorticity  may  be 
entrained  into  the  recirculation  cell,  thereby  reducing  the  total  circu¬ 
lation.  The  entrainment  and  detrainment  of  the  recirculation  cell  have 
been  discussed  in  detail  by  Crow  [Ref.  51]  and  by  Tombach  [Ref.  59], 
among  others,  and  will  not  be  discussed  here  further.  Suffice  it  to 
note  that  the  end  result  of  the  aging  of  vortices  is  to  decrease  their 
upward  migration. 

The  next  demise  mechanism  is  that  first  rigorously  analyzed  by 
Crow  [Ref.  54]  and  known  either  as  Crow  instability  or  sinusoidal  insta¬ 
bility.  Crow  has  shown  that  two  initially  parallel  trailing  vortices  are 
unstable  to  small  perturbations  and  links  at  intervals  of  8.6  times  the 
initial  core  spacing  bQ.  The  rings  so  formed  gradually  dissipate  at  a 
rate  dependent  on  atmospheric  conditions.  In  fact.  Crow  and  Bate  [Ref. 

55]  have  shown  the  role  played  by  turbulence  in  the  sinusoidal  instability, 
the  wavelength,  and  the  demise  of  the  entire  trailing  vortex.  In  spite  of 
the  uncertainties  associated  with  the  quantification  of  the  effects  of 
turbulence,  the  instability  analysis  of  Crow  still  remains  as  the  only 
rigorous  mathematical  model  of  any  one  of  the  demise  mechanisms. 

In  the  data  shown  in  Fig.  32,  the  Crow  instability  was  observed 
in  the  trailing  vortices  for  stratification  parameters  less  than  about 
0.75.  At  higher  stratification  parameters,  the  vortices  do  not  last  long 
enough  to  experience  Crow  instability.  This  raises  the  important  ques¬ 
tion  as  to  the  time  required  from  the  inception  to  the  linking  of  the 
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vortices  in  stratified  fluids.  Such  an  analysis  does  not  yet  exist  and 
the  said  time  will  probably  be  a  strong  function  of  the  ambient  turbu¬ 
lence  as  in  the  case  of  unstratified  fluids. 

The  events  following  the  linking  of  the  vortices  are  both  inter¬ 
esting  and  complex.  The  stretching  of  the  vortices,  the  instability 
propagating  about  the  rings,  the  increase  of  the  core  size,  and  the 
generation  of  opposite-signed  vorticity  give  rise  to  as  yet  incalculable 
consequences.  This  aspect  of  the  investigation,  like  many  other  inter¬ 
esting  problems  raised  during  the  course  of  this  investigation,  needs 
further  examination,  probably  through  the  use  of  high-speed  motion 
pictures. 

The  last  and  certainly  one  of  the  most  important  demise  mechanisms 
to  be  discussed  is  the  vortex  breakdown.  When  the  axial  velocity  in  a 
vortex  becomes  nearly  equal  to  the  maximum  tangential  velocity  at  the  edge 
of  the  core,  the  vortex  undergoes  a  rapid  transition  from  supercritical 
to  subcritical  flow.  This  abrupt  transition  may  occur  in  one  of  three 
forms  depending  on  the  circulation  and  the  Reynolds  number  of  the  ambient 
flow.  These  are  axisymmetric  vortex  breakdown,  spiral  breakdown,  and 
double-helix  breakdown  [Ref.  7].  The  occurrence  of  the  vortex  breakdown 
may  be  precipitated  by  an  adverse  pressure  gradient.  For  example,  the 
swirling  flow  in  a  conical  diffuser  produces  vortex  breakdown  at  a  dis¬ 
tance  must  closer  to  upstream  than  a  uniform  pipe.  Thus  in  the  case  of 
an  axisymmetric  streamlined  body,  t  ^  deceleration  of  the  ambient  flow  at 
the  aft  end  produces  vortex  breakdown  in  the  trailing  vortex.  In  the  case 
of  the  Delta-wings,  the  breakdown  results  from  the  expansion  of  the  core 
and  the  reduction  of  the  tangential  velocity.  The  consequence  of  the 
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breakdown  is  not  the  total  annihilation  of  the  vortex,  as  it  was  once 
thought.  Observations  and  measurements  have  shown  that  the  vortices  sur¬ 
vive  the  vortex  breakdown  and  continue  to  rise.  It  is  not  yet  known  as 
to  how  much  circulation  if  left  in  the  vortex  following  the  vortex  break¬ 
down  or  a  cascade  of  breakdowns.  In  all  cases,  the  trailing  vortex  may 
contain  one  or  more  velocity  maximums  with  numerous  smaller  swirls.  It 
appears  that  the  reduction  of  circulation  is  more  with  vortex  breakdown 
than  with  the  other  mechanisms. 

The  data  points  shown  in  Fig.  32,  particularly  for  small  values 
of  the  stratification  parameter  have  resulted  from  the  vortices  which  had 
sufficient  time  to  undergo  Crow  instability  and  vortex  breakdown.  Depend¬ 
ing  on  the  particular  experiment,  either  the  Crow  instability  or  the  vortex 
breakdown  occurred  first  since  there  was  no  strong  adverse  pressure  gradi¬ 
ent  to  precipitate  the  vortex  breakdown  first.  Since  the  vortex  motion 
resulting  from  a  combination  of  Crow  instability  followed  by  vortex  break¬ 
down  is  not  the  same  as  that  resulting  from  vortex  breakdown  followed  by 
some  form  of  Crow  instability,  one  does  expect  considerable  scatter  in 

H  ,.,/b  .  This  is  evidenced  by  the  data  points  in  Fig.  32,  particularly 
max  o 

for  Nbo/VQ  =0.  It  must  be  emphasized  that  Hmax/bQ  is  a  measure  of  the 
maximum  relative  height  at  which  the  vortex  cores  have  essentially  dissi¬ 
pated.  It  is  quite  possible  that  some  residual  vorticity  may  still  per¬ 
sist  in  the  flow.  This  vorticity  is  not  expected  to  have  significant 
consequences  relative  to  the  vorticity  associated  with  the  trailing 
vortices  prior  to  the  dissipation  of  the  cores.  In  other  words,  the 
vorticity  density  in  the  flow  field  for  H/b„  larger  than  H_av/b_  is  insigni- 
ficant  relative  to  H/bQ  smaller  than  Hmax/b0 • 
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The  need  to  determine  the  lifespan  of  vortices  and  to  understand 
the  physics  of  the  demise  mechanisms  gave  rise  to  several  approximate 
analyses.  These  will  now  be  discussed  in  some  detail  and  another 
approximate  analysis  will  be  presented. 

Scorer  and  Davenport  [Ref.  50]  assumed  that  the  impulse  imparted 
to  the  vortex  pair  and  the  circulation  of  each  vortex  remain  constant. 
This  resulted  in  a  gradually  decreasing  distance  between  the  vortices 
given  by: 


H  «  jgj-  sinh  (ANt) 


(8) 


and 


b 


b 


1 

o  cosh  (ANt) 


(9) 


in  which 


A  =  0.67 


Eliminating  time  between  Eqs.  (8)  and  (9),  one  has 


/,b o,2  ,  ANH 

V(v>  -  yr 


(10) 


Scorer  and  Davenport  assumed  that  the  recirculation  cell  will  cease  to 
rise  when  the  vortices  annihilate  each  other  by  having  their  cores  touch. 
In  other  words,  assuming  b  =  2rQ,  one  has  from  Eq.  (10), 
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The  above  equation  yields  the  maximum  relative  height  in  terms  of  the 
stratification  parameter  provided  that  one  has  some  idea  about  the  core 
radius,  or  b0/rQ.  The  use  of  rQ  =  0 - 1 bQ  in  Eq.  (11)  results  in  unrealis¬ 
tically  large  Hmax/bQ  values.  Various  values  of  bQ/r0  were  examined  so 
as  to  obtain  some  correspondence,  at  least  for  some  ranges  of  NbQ/V0. 

This  resulted  in  a  choice  of  bQ/r0  *  3.  Then  Eq.  (11)  reduces  to 
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Equation  (12)  is  plotted  in  Fig.  33  together  with  the  experimental  data. 
Clearly,  there  are  two  fundamental  difficulties  associated  with  Scorer 
and  Davenport's  analysis.  First,  the  assumed  core  radius  of  bQ/3  is  too 
large  to  be  realistic.  Second,  as  the  stratification  parameter  approaches 
zero,  the  analysis  predicts  increasingly  large  ,/b„  values,  implying 
that  a  trailing  vortex  in  an  unstratified  infinite  fluid  medium  will  never 
cease  to  rise.  These  shortcomings  are  simple  consequences  of  the  fact 
that  the  physics  associated  with  the  demise  mechanisms  have  not  been  in¬ 
corporated  into  the  analysis.  Judging  from  the  results  of  other  works 
since  the  pioneering  work  of  Scorer  and  Davenport,  one  must  state  in  all 
fairness  that  it  has  not  yet  been  possible  to  incorporate  the  consequences 
of  turbulence.  Crow  instability,  and  vortex  breakdown  into  the  analysis. 
Attempts  made  to  model  numerically  to  include  the  effect  of  turbulence 
will  be  discussed  later. 
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Experimental  D 


Another  approximate  analysis  of  the  ultimate  vertical  migration 
of  the  recirculation  cell  may  be  carried  out  by  assuming  that  the 
initial  kinetic  energy  of  the  fluid  within  the  cell  is  converted  to 
potential  energy. 

The  kinetic  energy  of  the  fluid  within  the  recirculation  cell 
is  given  by  [Ref.  62] 

E  =  ^pr2(ln  |)  (13) 

The  potential  energy  of  the  reel rculation  cell  may  be  written  as 

p  •  <14> 

2 

in  which  mbQ  is  the  area  of  the  recirculation  cell  and  m  =  2.83.  Equat¬ 
ing  the  kinetic  and  potential  energies  using  rQ/bQ  *  0.1,  and  simplifying 
one  has 

Hmax/bo  =  (Nb’/V0)  (15) 

Equation  (15)  is  shown  in  Fig.  33  together  with  Eq.  (12)  and  experi¬ 
mental  data.  Clearly,  Eq.  (15)  represents  the  upper  bound  of  the 
Hmax/^o  values  since  no  energy  dissipation  due  to  any  one  of  the  demise 
mechanisms  has  been  incorporated  into  the  analysis.  Finally,  it  should 
be  emphasized  that  the  predictions  of  Eq.  (15)  could  have  been  partially 
matched  to  the  experimental  data  by  artificially  increasing  the  relative 
core  radius  rQ/bQ.  In  the  foregoing,  r0/bQ  was  assumed  to  have  the 
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experimentally  observed  value  of  0.1  in  order  to  accentuate  the  energy 
loss  brought  about  by  the  complex  demise  mechanisms. 

A  third  approximate  analysis  may  be  carried  out  by  assuming  that 
the  vortex  separation  remains  constant,  but  the  strength  of  the  vortices 
vary  with  time. 

Starting  with  the  analysis  of  Scorer  and  Davenport  [Ref.  50], 
i.e.,  by  equating  the  time  rate  of  change  of  impulse  imparted  to  the 
vortex  pair  to  the  buoyant  force,  one  has 


3(prbQ) 

at 


g  m  bQ  Ap 


(16) 


Assuming  3bQ/3t  =  0,  one  has 


pbo  ar  3  9  m  bjj  p  B  z 

in  which  Ap  =  p B z  and  B  =  N  /g.  Noting  that  V  = 
may  be  reduced  to 


(17) 


r 

2lTb0 


,  Eq.  (17) 


=  gmb0Bzdz  (18) 

o 

Integration  of  Eq.  (18)  and  rearrangement  in  terms  of  the  stratification 
parameter  yields 


Hmax  _  1.49 


(19) 


The  correspondence  between  the  experimental  data  and  the  prediction  of 
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Eq.  (19)  (see  Fig.  34)  may  be  somewhat  fortuitous  because  of  a  number  of 
reasons.  The  size  of  the  vortex  core  does  not  enter  into  the  analysis 
even  though  it  is  known  to  have  an  important  effect  on  the  life  of 
vortices.  Furthermore,  as  in  all  other  analyses,  the  known  demise  mechan¬ 
isms  do  not  appear  in  Eq.  (19). 

All  of  the  foregoing  approximate  analyses  point  out  emphatically 
that  there  are  no  simple  approximations  with  which  one  can  predict  the 
rise  of  the  vortices.  One  must  use  experimental  data  and  heuristic 
reasoning  to  delineate  the  role  played  by  various  demise  mechanisms  in 
the  ascent  of  vortices.  This  fact  is  further  accentuated  by  Fig.  35 
where  the  experimental  data  are  shown  together  with  a  best  fit  curve 
defined  by 


4.64  exp[-1.3(NbQ/Vo)2] 


(20) 


The  particular  dependence  of  Hmax/b0  on  NbQ/Vo  strongly  suggests  that  the 
phenomena  such  as  turbulence,  entrainment  and  detrainment.  Crow  instability, 
and  the  vortex  breakdown  play  far  more  complex  roles  than  what  may  be 
quantified  in  a  simple  approximate  analysis. 


E.  SAILPLANE  GENERATED  VORTICES 

The  results  obtained  with  a  streamlined  body  and  its  lifting  surface 
in  stratified  and  unstratified  water  will  be  presented  in  the  following. 

Figure  36  shows  the  normalized  vertical  distance  H/bQ  as  a  function 
of  the  relative  distance  Z'/bQ  along  the  direction  of  motion  for  three 
representative  model  velocities.  The  reason  for  the  choice  of  the  particu¬ 
lar  normalizing  parameters  is  that  the  body  proximity  does  not  allow 
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Sailplane  Generated  Vortices  as  a  Function 


the  determination  of  an  initial  rise  velocity  VQ.  In  spite  of  the 
scatter,  the  data  show  that  there  is  very  little  or  no  dip  in  the  path 
of  the  trailing  vortices  as  they  approach  the  aft  end  of  the  body  at 
2VbQ  a  14.  Once  the  vortices  are  free  from  the  body  proximity  effects, 
they  continue  to  rise  at  a  nearly  constant  velocity.  In  fact  it  is  this 
particular  rise  velocity  that  is  used  in  Fig.  37  to  plot  H/bQ  as  a  function 
of  V0t/bQ. 

Also  shown  in  Fig.  36  is  the  analytical  result  based  on  ideal  line 
vortices  (see  Appendix  E).  The  vortex  characteristics  used  in  the  calcu¬ 
lations  were  those  appropriate  for  the  lifting  surface  of  the  streamlined 
body.  The  predictions  of  the  inviscid  flow  model  come  surprisingly  close 
to  the  experimental  data.  This  may  be  due  to  a  number  of  reasons,  the  most 
important  of  which  is  the  use  in  the  calculations  of  the  vortex  strength 
which  corresponded  to  that  experimentally  determined  value  after  a  series 
of  vortex  core  bursts,  i.e.,  the  use  of  the  vortex  velocity  based  on  the 
data  shown  in  Fig.  36. 

The  exploratory  results  obtained  with  one  stratification  parameter 
(SP  =  0.625)  are  shown  in  Figs.  38  and  39.  A  comparison  of  these  figures 
with  Figs.  36  and  37  show  that  the  stratification  inhibits  the  rise  of 
the  vortices  as  anticipated  on  the  basis  of  the  experiments  with  Delta- 
wings.  The  effects  of  other  stratifications  on  the  evolution  of  trailing 
vortices  shed  from  the  lifting  surfaces  of  an  axi symmetric  slender  body 
will  have  to  be  examined  in  light  of  additional  experiments.  The  determina¬ 
tion  of  the  circulation  existing  in  the  trailing  vortices  after  the  vortex 
breakdown  will  enable  one  to  compare  with  greater  confidence  the  prediction 
of  the  numerical  models  with  the  results  obtained  experimentally. 
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Figure  37,  The  Rise  of  the  Sailplane  Generated  Vortices  as  a  Function  of  V  t/b 
for  SP  =  0. 
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Figure  38.  The  Rise  of  the  Sailplane  Generated  Vortices  as  a  Function  of  Z/b( 
for  SP  =  0.625. 


F.  NUMERICAL  MODEL  OF  LAMINAR  VORTEX  PAIR 


As  discussed  in  the  previous  sections,  the  need  to  determine  the 
lifespan  of  vortices  and  to  understand  the  physics  of  the  demise  mechan¬ 
isms  gave  rise  to  several  approximate  analyses.  Numerical  methods  have 
been  employed  only  recently  using  computers  to  attempt  to  incorporate  the 
effect  of  turbulence  in  the  evolution  and  ultimate  demise  of  the  vortex 
pair.  It  is  hoped  that  such  methods  will  provide  not  only  a  power  of 
prediction  of  the  migration  of  trailing  vortices  but  also  some  insight  as 
to  the  relationship  between  turbulence  and  the  demise  mechanisms.  To  be 
sure,  there  are,  at  present,  no  numerical  methods  dealing  with  the  three- 
dimensional  trailing  vortices  and  the  effects  of  vortex  breakdown  and 
Crow  instability. 

The  most  advanced  numerical  methods  have  been  applied  to  the  study 
of  the  two-dimensional,  unsteady  vortex  pair.  Hecht  et  al.  [Ref.  44], 
introduced  such  a  model  which  solves  the  mean  and  ensemble-averaged 
Reynolds  stress  equations  of  fluid  motion,  based  on  a  second-order  closure 
turbulence  model.  The  initial  distribution  of  vorticity  was  assumed  to 
be  Gaussian  of  the  form 

C  *  exp[-(r/o)^]  (21) 

TO 

with  o  *  0.25bQ  (this  corresponds  to  rQ/b0  a  0.177).  The  initial  turbu¬ 
lence  distribution  was  again  Gaussian  with  a  spread  of  0.25bQ,  and  was 
assumed  to  be  Isotropic  with  a  representative  integral  scale.  An  upward 
velocity  was  employed  to  follow  the  vortex  descent.  A  Reynolds  number 

5 

of  V  b  /v  *  2*10  was  assumed  for  the  analysis, 
oo 


The  model  was  validated  by  comparing  the  calculated  and  experimentally 
determined  trajectories  of  vortex  rings.  The  results  of  this  model  for 
a  vortex  pair  in  atmosphere  with  a  stratification  parameter  of  0.8  are 
compared  In  Fig.  40  with  those  of  the  Delta-wings  In  water  with  a  strati¬ 
fication  parameter  of  0.75.  The  model  matched  the  experimental  data  very 
well,  although  there  are  no  physical  arguments  to  justify  the  selected 
integral  scale  for  turbulence.  Furthermore,  the  value  of  a  should  be 
about  half  that  assumed. 

In  a  recent  paper,  Hecht,  Bilanin,  and  Hirsh  [Ref.  53]  compared  the 
predictions  of  the  aforementioned  model  with  other  analytical  [Refs.  49, 

51,  and  59]  and  experimental  results  [Ref.  43].  Initial  distribution  of 
vorticity  was  based  on  the  data  from  Burnham  et  al.  [Ref.  43].  The  initial 
a  used  in  Eq.  (21)  was  determined  to  be  0.074bQ.  The  integral  scale  for 
turbulence  was  set  to  a  smaller  value  for  this  analysis  and  a  sensitivity 
study  of  the  macroscale  effects  on  vortex  descent  was  presented.  The 
results  of  this  model  for  SP  =  0.67  and  a/bQ  *  0.075  are  compared  with  the 
Delta-wing  results  for  SP  =  0.625  and  SP  -  0.75  in  Fig.  41.  This  model 
does  not  match  the  data  for  VQt/b0  greater  than  about  2. 

A  better  match  between  the  measured  and  calculated  values  could  have 

been  achieved  with  a  =  0.25bQ.  But  this  value  is  larger  than  that  in 

trailing  vortices  by  a  factor  of  about  3  relative  to  the  observations  of 

Burnham  et  al.  [Ref.  43].  Of  course  the  integral  scale  of  turbulence 

could  have  been  adjusted  also  to  achieve  more  realistic  results  for 

V„t/b„  larger  than  about  1.5. 
oo 

While  the  numerical  models  cited  above  have  pioneered  in  the  study  of 
turbulent  effects  in  trailing  vortices,  the  difficulties  associated  with 
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Comparison  ulth  Hecht  et  al  Model  C44D 


Experimental  Data  with  that  Predicted  by  Hecht  et  al  (Ref.  44) 


Figure  41,  Comparison  of  Experimental  Data  with  that  Predicted  by  Hecht  et  al  (Ref,  53), 


the  choice  of  free  parameters  (o  and  the  integral  scale  of  turbulence) 
point  out  the  need  for  additional  work. 

In  an  effort  to  eventually  produce  such  a  fully  turbulent  flow  model 
to  compare  the  three-dimensional  trailing  vortex  data,  gathered  from  the 
experimental  models  in  the  water  basin,  a  two-dimensional,  laminar,  un¬ 
steady  flow  numerical  model  of  the  vortex  pair  was  developed  as  part  of 
this  investigation. 

For  this  purpose,  the  Navier-Stokes  equations  were  written  as 


3M.  +  u  3*t  +  v 

at  u  ax  v  ay 


0-lf^vA 


+  3v  +  av 

at  u  ax  ay 


l  ap  .  .2 

9-  P  ay+vA  v 


(22) 

(23) 


where  the  y-axis  is  assumed  to  be  directed  downwards.  Eliminating  the 
pressure  between  the  two  equations,  one  has 


and  represents  the  vorticity.  For  flows  in  which  the  gravitational 
acceleration  is  several  orders  of  magnitude  larger  than  the  fluid  acceler¬ 
ations,  the  terms  in  the  bracket  in  Eq.  (24)  may  be  neglected.  This 
amount  to  invoking  the  Boussinesq  hypothesis.  Thus,  the  equation  of 
motion  and  the  equation  of  continuity  may  be  written  as 
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(26) 


,  u72p  (27) 

The  last  term  in  Eq.  (26)  represents  the  effect  of  the  density  gradient 
and  gives  rise  to  opposite-signed  vorticity. 

The  equations  that  were  chosen  to  describe  the  motion  include  a 
Gaussian  distribution  of  the  initial  vorticity  given  by 


C  *  iJ?exp[-(r2/2r2)]  (28) 

2-rr  rZ 
o 

Boundary  values  of  the  u  and  v  velocity  components  are  determined  by 
the  Biot-Savart  Integrals 


u(x,y) 


v(x,y) 


(29) 


(30) 


Where: 

r2  =  (x-x1)2  +  (y-y1)2 

x',y‘  denote  the  position  of  vorticity,  and 

x,y  the  point  where  the  velocity  components  are  calculated. 

The  stream- function  values  around  the  boundaries  and  velocity  within  the 
field  boundaries  are  calculated  from  u  and  v  using 
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u 


3£  . 

3y  * 


v 


9tj/ 

3x 


(31) 


Field  stream-function  values  are  calculated  using  Poisson's  equation 

v2<j,  =  -  Z  (32) 

The  calculations  had  to  be  carried  out  in  a  finite  numerical  mesh,  as 
in  all  other  numerical  analyses.  The  appropriate  boundary  conditions  are 
shown  on  the  boundaries  of  the  quadrant  in  which  the  calculations  have 
been  performed  (see  Fig.  42).  It  should  be  noted  that  the  conditions 
u(0,y)  =  v(x,0)  =  0  are  satisfied  automatically  by  evaluating  the  Biot- 
Savart  equations  in  all  four  quadrants.  Furthermore,  the  symmetry  pro¬ 
vided  by  the  free  surface  (x-axIs)  and  by  the  normal  bisecting  the  vortex 
pair  (y-axis)  enable  one  to  confine  the  domain  of  calculation  to  a 
single  quadrant. 

Additional  details  of  the  analysis  are  as  follows: 

(i)  The  domain  of  calculation  is  prescribed; 

(ii)  The  vorticity  distribution  given  by  Eq.  (28)  is  assigned  to 
mesh  points; 

(iii)  A  particular  stratification  is  assigned  to  the  density  distribution 

(iv)  Boundary  values  of  velocities  are  calculated  using  the  Biot- 
Savart  equations; 

(v)  The  boundary  values  of  the  stream- function  are  calculated; 

(vi)  The  flow  field  values  of  the  stream-function  are  calculated 
using  Poisson's  equation  and  a  successive  overrelaxation  scheme; 

(vii)  The  velocity  field  is  calculated  through  the  use  of  the  stream- 
function  values. 
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Figure  42.  Computational  Field  and  Boundary  Conditions  for 
the  Numerical  Model . 
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The  foregoing  marks  the  end  of  the  first  time  step.  The  successive 
steps  are  sequenced  as  follows: 

(i)  The  vorticity  and  density  fields  are  calculated  using  an 
upwind-differencing  scheme; 

(ii)  The  new  center  of  vorticity  is  calculated  using  the  square  of 
vorticity  as  the  weighing  parameter; 

(iii)  A  vertical  velocity  is  imposed  on  the  field  using  the  velocity 
calculated  from  the  rate  of  change  of  the  center  of  vorticity; 
and  finally, 

(iv)  The  new  points  at  the  top  of  the  field  are  initialized  and  the 
vorticity  is  used  to  calculate  the  boundary  velocities. 

The  last  step  completes  the  iteration  and  marks  the  beginning  of  a  new 
sequence. 

The  plots  of  the  streamlines  of  the  flow  field  for  a  simulation  of 

the  motion  of  the  vortices  generated  by  Delta  1  for  SP  =  0  are  shown  in 

Figs.  43  through  47.  The  time  indicated  on  each  plot  is  the  normalized 

time  t/(bQ/Vo)  which  is  of  course  identical  to  its  common  form  v0t/bQ. 

In  these  calculations,  D/bQ  was  taken  to  be  4.31,  and  the  initial  value 

of  r  was  assumed  to  be  0.09b  .  The  calculated  values  of  H/b„  are  com- 
oo  o 

pared  with  those  obtained  experimentally  in  Fig.  48.  The  matching  of 
the  two  results  is  surprising  in  view  of  the  fact  that  neither  turbulence 
nor  any  other  demise  mechanism  in  conjunction  with  turbulence  was  taken 
into  consideration.  The  surprising  agreement  between  the  two  results 
may  be  attributed  to  numerical  diffusion  and  to  the  relative  stability 
of  the  trailing  vortices  in  unstratified  fluids. 

Sample  streamline  plots  resulting  from  the  numerical  calculation  for 
SP  -  0.50  are  shown  in  Figs.  49  through  53  for  various  values  of  H/bQ. 
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Figure  43.  Streamlines  of  the  Numerical  Model  for  SP 
H/b0  =  0. 
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Figure  44.  Streamlines  of  the  Numerical  Model  for 


SP  =  0,  H/b0  =  0.60. 
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Figure  45.  Streamlines  of  the  Numerical  Model  for 
SP  =  0,  H/bQ  =  1.13. 
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Figure  46.  Streamlines  of  the  Numerical  Model  for 
SP  =  0,  H/b  =  1.59. 
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Figure  47.  Streamlines  of  the  Numerical  Model  for 
SP  =  0,  H/bQ  =  2.01. 
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Figure  48.  Comparison  of  Experimental  Data  with  the  Predictions  of  the  Numerical 
Model  for  SP  =  0. 
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Figure  49.  Streamlines  of  the  Numerical  Model  for 
SP  =  0.50,  H/b  =  0. 
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Figure  50.  Streamlines  of  the  Numerical  Model  for 
SP  =  0.50,  H/bQ  =  0.61. 
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Figure  51.  Streamlines  of  the  Numerical  Model  for 
SP  =  0.50,  H/bQ  =  1.17. 
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Figure  52.  Streamlines  of  the  Numerical  Model  for 


SP  =  0.50,  H/bQ  =  1.70. 
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Figure  53.  Streamlines  of  the  Numerical  Model  for 
SP  =  0.50,  H/bQ  =  2.22. 
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The  comparison  with  Figs.  43  through  47  show  that  in  the  stratified  case 
the  vortices  do  not  split  out  laterally  as  much  as  in  the  unstratified 
case.  The  reason  for  this  is  that  the  opposite-signed  vorticity  is 
generated,  as  shown  by  numerical  calculations,  outside  the  lower  right 
and  lower  left  sides  of  the  recirculation  cell.  The  effect  of  this 
vorticity  is  to  pull  the  vortices  together  and  thus  increase  the  velocity 
of  their  upward  migration.  Had  turbulent  diffusion  been  present,  the 
reduction  in  circulation  would  have  more  than  offset  the  effect  of  opposite- 
signed  vorticity,  thereby  slowing  the  rise  of  vortices.  In  laminar 
vortices,  however,  the  strength  of  the  vortices  changes  very  slowly  and 
the  vortices  accelerate  upward  in  an  amount  closely  related  to  the  magni¬ 
tude  of  stratification.  The  foregoing  arguments,  however  qualitative, 
point  out  once  again  the  conflicting  conclusions  arrived  at  by  various 
workers  (see  Table  I) . 

The  experimental  and  calculated  values  are  compared  in  Figs.  54  and 
55  for  SP  =  0.50  and  SP  =  1.0,  respectively.  The  correspondence  between 
the  measured  and  calculated  values  is  surprisingly  good,  partly  because 
of  the  fact  that  even  the  experimental  values  of  H/bQ  for  the  stratified 
(SP  =  0.50)  and  unstratified  cases  do  not  significantly  differ  (see  Figs. 

23  and  27).  However,  in  the  case  of  SP  =  1.0  where  the  experimental  data 
for  SP  =  0  and  SP  =  1 .0  show  significant  differences  (see  Figs.  23  and  27), 
the  predictions  of  the  numerical  calculation  differ  significantly  from 
the  experimental  data  as  seen  from  Fig.  55.  As  noted  earlier,  the  strong 
opposite- signed  vorticity  generated  in  this  case  accelerates  the  vortices 
upward  and  yields  larger  H/bQ  values.  It  is  interesting  to  note  that  at 
the  early  stages  of  the  motion,  where  the  effect  of  counter-signed 
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Figure  54*  Comparison  of  Experimental  Data  with  the  Predictions  of  the  Numerical 
Model  for  SP  =  0.50. 
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Figure  55.  Comparison  of  Experimental  Data  with  the  Predictions  of  the  Numerical 
Model  for  SP  =  1.0. 


vorticity  is  not  yet  evident,  the  slope  of  the  numerical  prediction  is 
nearly  identical  to  that  of  the  experimental  data. 

It  should  be  emphasized  that  the  foregoing  was  only  an  exploratory 
numerical  investigation  of  the  evolution  of  a  laminar  vortex  pair.  It 
served  to  confirm  the  need  for  the  development  of  a  comprehensive 
numerical  model  which  is  not  only  true  to  the  underlying  physics  of  the 
phenomena  but  also  devoid  of  floating  parameters  whose  values  are  arbitrarily 
chosen  (e.g.,  initial  vorticity  diffusion,  rQ/b0 ,  mesh  size,  size  of 
calculation  domain,  etc.).  The  results  also  show  that  turbulence  plays 
an  extremely  important  role  in  the  diffusion  of  vorticity  and  thus  in 
the  entire  life  of  vortices  from  their  inception  to  their  ultimate  demise. 

Even  when  such  a  comprehensive  model  is  developed  it  does  not  now 
appear  that  it  would  be  too  easy  to  surmount  the  difficulties  associated 
with  the  modeling  of  the  consequences  of  the  Crow  instability,  vortex 
breakdown,  and  the  entrainment  and  detrainment  of  the  recirculation  cell. 

Of  course  it  is  not  too  surprising  that  such  an  investigation  should 
require  a  better  understanding  of  turbulence  itself  through  measurements 
and  calculations  of  vortical  flows. 

G,  SCALE  EFFECTS 

The  results  obtained  with  Delta  wings  have  shown  that  the  Reynolds 
number  does  not  play  a  measurable  role  in  the  evolution  of  the  turbulent 
trailing  vortices.  Thus,  the  data  presented  herein  may  be  used  for  larger 
Delta  wings  provided  that  the  background  turbulence  Is  negligible  and  the 
stratification  is  linear.  Additional  experiments  are  needed  to  delineate 
the  effects  of  the  background  turbulence  and  the  nonlinearity  of  the 
stratification.  Such  experiments  are  currently  underway. 
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As  to  the  streamlined  body,  the  Reynolds  number  of  the  model  is  about 
800  times  smaller  than  that  of  a  full  scale  vehicle.  In  spite  of  the 
rather  large  difference  between  the  Reynolds  numbers  of  the  model  and  the 
prototype,  however,  the  data  presented  herein  are  thought  to  be  valid  for 
the  full  scale  vehicle  provided  that  the  background  turbulence  is  negli¬ 
gible  and  the  stratification  is  nearly  linear.  This  assertion  is  based 
on  the  fact  that  the  boundary  layers  on  the  model  and  its  lifting  surfaces 
were  rendered  turbulent  with  tripping  wires  and  sand  strips.  Obviously, 
this  does  not  preclude  the  need  for  full  scale  experiments  for  a  number 
of  reasons,  the  most  important  ones  being  the  effect  of  background  turbu¬ 
lence  and  the  nonlinearity  of  the  stratification. 

One  must  bare  in  mind  that  there  are  fundamental  differences  between 
the  demise  mechanisms  of  trailing  vortices  generated  by  a  Delta  wing  or 
rectangular  wing  (in  the  absence  of  any  body-proximity  effects)  and  those 
generated  by  a  vehicle  and  its  control  surfaces.  The  demise  of  the  vortices 
emanating  from  a  Delta  wing  is  primarily  due  to  sinusoidal  instability  and 
turbulent  diffusion.  Thus,  the  lifespan  of  these  vortices  may  be  strongly 
affected  by  the  background  turbulence  [Ref.  55].  The  demise  of  the  trailing 
vortices  generated  by  the  lifting  surfaces  of  an  axisymmetric  body  is 
primarily  due  to  the  vortex  breakdown.  In  fact,  the  data  presented  herein 
show  that  the  vortex  breakdown  (brought  about  by  the  body-proximity  effects) 
severely  limits  the  rise  of  the  trailing  vortices.  The  comparison  of  the 
data  obtained  with  Delta  wings  and  the  streamlined  body  show  that  the  maximum 
rise  of  the  cores  of  the  vortices  generated  by  the  sail  planes  of  the  stream¬ 
lined  body  is  about  2.5  times  smaller  than  that  of  the  Delta  wings.  Thus, 
in  the  case  of  vehicles  one  would  be  interested  not  only  with  the  effect  of 
background  turbulence  and  the  type  of  stratification  but  also  with  the 
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effect  of  Reynolds  number  on  the  intensity  and  generation  of  vortex 
breakdown  or  cascade  of  vortex  breakdowns.  The  previous  studies  [Refs. 

6-8]  have  shown  that  the  occurrence  of  vortex  breakdown  is  not  materially 
affected  by  the  Reynolds  number  (and  Mach  number)  provided  that  the 
vortex  is  initially  turbulent.  In  other  words,  the  background  turbulence 
is  not  of  major  importance  in  the  occurrence  of  vortex  breakdown,  at 
least  when  the  breakdown  takes  place  on  a  single  vortex.  One  cannot 
ignore  the  fact  that  events  such  as  vortex  breakdown  do  not  occur  simul¬ 
taneously  and  symmetrically  on  a  pair  of  trailing  vortices.  The  mutual 
induction  of  vortices  and  the  non-symmetrical  occurrence  of  cascades  of 
vortex  breakdown  on  each  vortex  may  enhance  the  role  played  by  the  back¬ 
ground  turbulence.  Furthermore,  the  amount  of  circulation  remaining  in 
the  recirculation  cell  following  the  breakdown  may  be  affected  by  the 
background  turbulence.  The  best  estimate,  at  present,  is  that  the  inten¬ 
sity  of  the  adverse  pressure  gradient  (i.e.,  the  mean  flow  field  about  the 
body)  is  fundamentally  responsible  for  the  occurrence  and  the  position  of 
the  vortex  breakdown.  The  subsequent  demise  of  the  vortices  may  be  affected 
by  the  background  turbulence,  mutual  induction,  and  the  type  of  stratifi¬ 
cation. 

It  is  evident  from  the  foregoing  that  it  would  be  very  desriable  to 
conduct  a  series  of  field  experiments.  In  doing  so,  the  data  presented 
herein  should  serve  well  for  their  planning  and  for  the  interpretation  of 
the  environmental  effects  on  the  evolution  of  the  trailing  vortices  and 
vortex  rings.  Clearly,  meaningful  full  scale  experiments  require,  for 
their  interpretation,  not  only  the  measurement  of  the  characteristics  of 
vortices  but  also  the  quantification  of  the  prevailing  environmental 
conditions. 
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In  summary,  the  evolution  of  vortices  shed  by  a  body  and  its  control 
surfaces  in  a  natural  environment  with  non-isotropic  background  turbulence, 
currents,  shear,  wave  motion,  and  stratification  is  such  a  complex  time- 
dependent  phenomenon  which  demands  quantitative  analysis.  Ingenious 
experiments,  and  qualitative  descriptions  for  the  purpose  of  explaining 
the  underlying  causes  of  the  physical  effects  and  providing  a  model  by 
which  design  and  operation  can  be  safely  and  reasonably  performed. 
Ultimately,  the  flow  field  and  the  vortex  motion  about  full  scale  vehicles 
will  not  be  quantified  solely  due  to  the  knowledge  gained  through  the  study 
of  well -posed  problems  and  controlled  laboratory  experiments  alone  but 
rather  through  the  interaction  and  reasoble  marriage  of  the  two  approaches: 
basic  studies  and  full  scale  experiments. 
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V.  CONCLUSIONS 


The  investigation  reported  herein  warranted  the  following  conclusions: 

1.  In  unstratified  water,  the  trailing  vortices  generated  by  a 
lifting  surface  moving  at  a  negative  angle  of  attack  rise  with  time  and 
dissipate  ultimately  at  a  finite  height  (about  H/bQ  *  4.6).  The  top  of 
the  recirculation  cell  rises  to  a  normalized  height  of  about  H/bQ  -  6. 

This  height  should  not  be  confused  with  the  height  to  which  the  dye  in 
the  trailing  vortex  diffuses  after  a  long  period  of  time. 

2.  The  rate  of  rise  of  vortices  is  always  smaller  than  that  predicted 
theoretically  through  the  use  of  the  inviscid  line  vortices. 

3.  The  rise  of  the  vortices  was  shown  to  be  governed  by  the  dimension¬ 
less  parameters  VQt/b0,  D/b0,  and  Nb0/VQ  (the  stratification  parameter)  for 
linearly  stratified  fluids  with  no  background  turbulence. 

4.  The  effect  of  stratification  is  to  reduce  the  rate  of  rise  of 
vortices  and  the  maximum  height  attained  by  them.  In  fact,  the  larger  the 
stratification  (i.e.,  NbQ/V0),  the  smaller  the  maximum  height  attained  by 
the  vortices. 

5.  The  ultimate  demise  of  the  vortices  generated  by  Delta  wings  is 
observed  to  be  brought  about  by  sinusoidal  instability  and  turbulent 
diffusion.  The  demise  of  the  vortices  generated  by  the  lifting  surfaces 
of  an  axi symmetric  body  is  brought  about  by  a  cascade  of  vortex  breakdowns 
and  the  turbulent  diffusion. 

6.  The  vortices  subjected  to  vortex  breakdown  rise  to  a  maximum 
height  which  is  about  2.5  times  smaller  than  those  subjected  to  sinusoidal 
instability. 

7.  The  results  have  shown  that  the  body-proximity  effects  on  trailing 
vortices  are  very  strong  and  always  give  rise  to  vortex  breakdown. 
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8.  Various  approximate  analytical  models  were  considered  and  all 
were  deemed  to  be  inadequate.  None  of  the  approximate  models  yielded  a 
finite  rise  height  in  unstratified  water  and  none  considered  the  effects 
of  the  demise  mechanisms. 

9.  An  exploratory  numerical  analysis  based  on  two-dimensional,  laminar 
unsteady  flow  assumptions  was  carried  out.  The  predictions  of  the  model 
were  in  good  agreement  in  the  unstratified  case.  In  the  stratified  case, 
the  generation  of  the  counter-signed  vorticlty  decreased  the  vortex  span 
and  increased  the  velocity  of  the  upward  migration  of  the  vortices.  The 
predictions  of  the  model  could  be  made  more  realistic  through  the  inclusion 
of  at  least  the  effects  of  turbulence  and  possibly  those  of  other  demise 
mechanisms. 

10.  The  scale  effects  were  examined  in  as  much  detail  as  possible. 

It  is  concluded  that  the  scale  effects  under  similar  conditions  (linear 
stratification  and  no  background  turbulence)  are  insignificant.  However, 
the  demise  of  the  vortices  undergoing  sinusoidal  instability  may  be  strongly 
affected  by  the  background  turbulence.  The  body-proximity  effects  and  the 
adverse  pressure  gradient,  rather  than  the  background  turbulence,  are  thought 
to  govern  the  occurrence  and  the  position  of  the  vortex  breakdowns.  However, 
the  ultimate  demise  of  the  vortices  subjected  to  vortex  breakdown  may  be 
affected  by  the  background  turbulence,  in  addition  to  other  laboratory  and 
environmental  factors  such  as  non-linear  stratification,  currents,  waves, 
wave-current  Interactions,  wakes  of  other  bodies,  shear,  etc. 
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APPENDIX  A 


RELATIONSHIP  BETWEEN  THE  AMOUNT  OF  SALT 
AND  STRATIFICATION  PARAMETER 


The  stratification  parameter  is  defined  as  NbQ/V0  and  is  given  by 


(A-l ) 


where: 

N  =  the  Brunt-Vaisala  frequency  and  is  defined  l>y 

«  •  (-^^)’/2  <A-2) 

bQ  =  the  vortex  core  separation 
2 

g  =  32.174  ft/s  is  the  gravitational  acceleration 

d  =  the  layer  thickness  in  inches 

P  =  the  fresh  water  density 
o 

AP  =  the  change  in  density  between  two  layers 
VQ  =  the  initial  vortex  rise  velocity. 


The  ratio  of  the  initial  vortex-rise  velocity  to  the  model  velocity  was 
found  to  be  constant  for  each  model,  independent  of  model  velocity,  i.e., 
VQ/U  =  Constant  for  a  given  model. 

The  brine-water  density  pp  in  the  n-th  layer  is  related  to  the 
concentration  Cn  by 
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p0(l  *  0.71C„) 


(A-3) 


where  the  salt  concentration  C  is  defined  in  general  by  the  following 
equation  as  a  function  of  the  amount  of  salt  added: 


1 

C 


u 

,  Fr. water 

+  w 

wsalt 


Thus,  one  has 


(A-4) 


AP 


°'71  WW 


or 


ap  - 

p 

o 


0.71 (Cn  - 


Cn-1  ^ 


( A-  5 ) 


and  the  Eq.  (A-l)  becomes, 

NbQ  bQ  J0.71xgx(Cn-Cj 

V0  "  (V0/U)U  >  d 


(A- 6) 


Solving  Eq.  (A-6)  with  respect  to  { Cn  -  C  ^ )  one  has 


C  -C  , 
n  n-1 


l077Tg 


(VU)U  p  Nb_  ; 

•)  -  [-4 — T  *  (-ir2)' 


(A- 7) 


where  n  =  1 ,2,3, . . .  . 

Because  of  the  fact  that  the  first  layer  normally  is  fresh  water 
without  salt,  i.e.,  CQ  =  0,  one  has 
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where  d  is  as  previously  described.  Therefore 


w 


n.salt 


561.6  x  d 


T 


nd  lr(»0/U)U  ,  Nb 


i  na  v  r ' 

T7^5T7T~g) 


J2(^)2 

o 


-  1 


For  example,  for 


NVVo  "  d  =  1"»  U  =  2.05  ft/sec,  and  VQ/U  * 

Wl,salt  3  0-01519  lbs  of  salt. 

Also,  for 


NVVo  *  1'°»  d  =  1"*  0  =2.05  ft/sec,  and  V  /U 

o 

W1  .salt  =  O-0607  lbs  of  salt,  and 
W2,salt  =  0.1215  lbs  of  salt. 


(A- 11 ) 


0.0391 


*  0.0391 
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APPENDIX  B 


TEMPERATURE  STRATIFICATION  EFFECTS 

A  nonuniform  temperature  stratification  which  might  result  from  the 
temperature  difference  between  the  ambient  air  and  the  water  in  the  basin 
could  have  an  adverse  effect  on  the  salt  stratification  and  on  the 
accuracy  of  the  data  obtained.  To  this  end,  an  analysis  has  been  carried 
out  by  assuming  that  the  basin  is  full  with  water  at  50  degree  temperature, 
the  basin  fills  at  a  much  faster  rate  than  the  rate  of  heat  diffusion, 
there  is  no  heat  transfer  between  the  bottom  and  side  walls  of  the  basin, 
and  the  air  temperature  at  the  air-water  interface  is  65  degrees. 

The  governing  equation  is  the  one-dimensional  unsteady  conduction 
equation,  given  by 


31  _ 

3t 


32T 


The  initial  and  the  boundary  conditions  are 


T(x ,0)  =  50°  F  =  T0 
T(L,t)  =  65°  F  =  T^ 

=  o 

3X 


Letting  9  =  T  -  Tw,  one  has 


39 

at 


329 

^7 


(B-l ) 


(B-2) 

(B-3) 

(B-4) 


(B-4) 
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and 


9(x,0)  *  e0  *  -  15°  F 

(B-5) 

9(L,t)  =  0 

(B-6) 

L9.(Q,t)  .  o 

3x 

(B-7) 

A  straightforward  application  of  the  technique  of  separation  of 
variables  yields 


T(x,t) 


T  +  2  9. 


00 


(-Dn 

V 


e 


t 

COS  X  X 

n 


(B-8) 


where 


,  _  (2n+l)ir 

An  2L 

For  the  case  under  consideration,  one  has 

t  =  4  hrs,  L  =  45  inches,  and  = 


(B-9) 


0.00022  in2/sec 


The  numerical  results  have  shown  that  the  temperature  distribution  in 
the  basin  is  affected  to  a  depth  of  only  seven  inches  after  a  period  of 
4  hours.  The  small  temperature  gradient  between  the  said  depth  and  the 
free  surface  gives  rise  to  a  negligible  density  gradient  in  all  but  the 
regions  very  close  to  the  free  surface,  a  region  which  is  not  of  special 
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interest.  It  is  of  importance  to  note  that  the  numerical  values  used 
in  the  calculations  correspond  to  those  encountered  in  coldest  months 
of  the  year.  For  the  majority  of  the  test  runs  the  water  temperature 
was  about  63  degrees  compared  to  the  air  temperature  of  67  degrees.  Thus, 
the  effect  of  heat  conduction  on  density  stratification  was  several  orders 
of  magnitude  smaller  than  that  calculated. 

As  to  the  effect  of  side  walls,  a  lateral  temperature  gradient  may 
develop  along  the  side  wall  since  the  thermal  conductivity  of  aluminum 
is  much  larger  than  that  of  water.  However,  such  a  temperature  gradient 
will  have  very  little  or  no  effect  on  the  vertical  gradient  and  hence  on 
the  motion  of  trailing  vortices.  Experiments  with  dye  particles  falling 
slowly  in  a  vertical  plane,  extending  from  wall  to  wall,  have  not  re¬ 
vealed  any  natural  convection  along  the  side  walls.  Thus,  it  is  concluded 
that  the  heat  transfer  in  any  direction  in  the  basin  did  not  play  an 
important'enough  role  to  affect  the  results. 


137 


APPENDIX  C  -  DELTA  1  TABULATED  DATA 
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\ 
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i 
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j 
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2.87 

3.46 

2.93 
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0.00 

0.00 
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0.00 
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.83 

.63  .34 

1.72 

1.36 
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.87  .78 
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1.93 
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1.88 
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3.26 
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1.81 
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1.78 
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4.32 

3.34 
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2.01 
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3.18 

2.14 
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£ 

1 

s 
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1.31 
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APPENDIX  E 


INTERACTION  OF  TIP  AND  ROOT  VORTICES  FOR 
VARIOUS  BODY  CONFIGURATIONS 


It  is  a  well-known  fact  that  the  combination  of  the  body  and  lifting 
surfaces  may  give  rise  to  root  vortices  in  addition  to  the  tip  vortices. 
The  occurrence  of  the  tip  vortices  depend  strongly  on  the  manner  the 
lifting  surfaces  are  attached  to  the  body.  If  and  when  they  occur,  they 
affect  the  evolution  of  the  tip  vortices  and  hence  their  ultimate  demise. 
Consequently,  it  is  of  special  importance  to  examine  in  detail,  albeit 
through  the  use  of  the  inviscid  line  vortices,  the  mutual  interaction  of 
the  two  types  of  vortices  for  various  body  shapes. 

The  analysis  described  below  assumes  three  different  body  shapes: 

(i)  a  streamlined  axisymmetric  body  with  lifting  surfaces;  (ii)  a  stream¬ 
lined  half  body  [half  of  that  used  in  (i)]  and  the  identical  lifting 
surfaces;  and  (iii)  only  the  lifting  surfaces  protruding  above  a  plane 
surface  ( representing  a  body  of  infinite  radius). 

The  case  (i)  was  analyzed  through  the  use  of  image  line  vortices  and 
the  appropriate  complex  velocity  potential  (see  Fig.  56).  The  analysis 
yields  the  following  complex  velocity 
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Figure  56.  Real  and  Image  Vortices  for  the 
Whole  Body  Analysis. 
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where  u  and  v  are  the  x  and  y  components  of  velocity;  and  the 
strengths  of  the  tip  and  root  vortices;  z ^  and  the  complex  displacement 
vector  locations  of  the  vortices  from  the  body  center;  a,  the  radius  of 
the  body  (dependent  on  the  distance  along  the  body);  and  A',  the  change 
in  body  cross-sectional  area  per  unit  longitudinal  length  of  the  stream¬ 
line  cylindrical  body. 

Writi  ng 
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where  At  is  the  time  step,  the  vortex  positions  were  calculated  as  a 
function  of  time.  Clearly,  the  values  superscripted  by  (n-1)  represent 
those  previous  to  those  superscripted  by  n.  The  value  of  A'  depends  on 
the  specific  shape  of  the  body  under  consideration. 

The  results  are  shown  in  Figs.  57  and  58  for  the  tip  and  root 
vortices.  In  these  figures,  Y/bQ  represents  the  normalized  height  of 
the  vortices  as  measured  from  the  streamlined  body  and  Z/bQ,  the  normalized 
distance  along  the  body,  starting  at  the  tip  of  the  lifting  surface. 

Figure  57  shows  the  vortex  positions  shortly  after  their  formation  while 
they  are  under  the  significant  influence  of  the  potential  flow  generated 
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Figure  57.  Positions  of  Tip  and  Root  Vortices  for  the  Whole  Body 
(Near  the  Lifting  Surface). 
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Figure  58.  Positions  of  Tip  and  Root  Vortices  for  the  Whole  Body 
(Away  from  the  Lifting  Surface). 


by  the  body.  It  should  be  noted  that  the  tip  vortex  dips  near  Z/bQ  3  14 
where  the  streamlined  body  terminates.  In  other  words,  the  potential 
flow  about  the  aft  end  of  the  body  pulls  the  tip  vortices  toward  the 
body.  Beyond  the  end  of  the  body,  the  four  vortices  are  under  each 
other's  mutual  induction  only,  without  body  effects.  Figure  58  shows 
the  long  term  interaction  between  the  tip  and  root  vortices  (of  course 
assuming  that  the  vortices  retain  their  inviscid  character).  Clearly, 
the  tip  and  root  vortices  rotate  about  each  other.  Real  vortices  are 
not  expected  to  survive  such  a  continuous  rotation. 

The  case  (ii),  (the  half-body  problem),  was  analyzed  in  the  same 
manner  through  the  use  of  the  vortex  system  shown  in  Fig.  59.  The 
following  complex  velocity  was  used  to  obtain  the  velocity  of  the  tip 
vortex  in  the  first  quadrant 
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in  which  every  parameter  has  the  previously  ascribed  meaning  to  it. 

The  results  of  the  half  body  analysis  are  presented  in  Figs.  60  and 
61.  The  vortex  positions  are  very  close  to  those  calculated  for  the 
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HALF  BODY  RNRLYSIS 
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Figure  60.  Positions  of  Tip  and  Root  Vortices  for  the  Half  Body 
(Near  the  Lifting  Surface). 


Figure  61.  Positions  of  Tip  and  Hoot  Vortices  for  the  Half  Body 
(Away  from  the  Lifting  Surface). 


whole  body,  particularly  shortly  after  the  formation  of  the  vortices. 

Thus,  one  may  conclude  from  a  hydrodynamic  as  well  as  practical  point 
of  view  that  the  half  body  closely  models  the  whole  body  while  providing 
additional  experimental  flexibilities  (e.g.,  larger  depth  to  body  diameter, 
mounting  of  the  test  body,  etc.). 

The  plane-wall  analysis,  i.e.,  case  (i i i ) ,  was  carried  out  through  the 
use  of  simple  images  as  shown  in  Fig.  62.  The  following  complex  velocity 
describes  the  velocity  of  the  tip  vortex  in  the  first  quadrant: 
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Once  again,  identical  values  of  the  parameters  were  used  in  order  to 
compare  the  results  of  the  three  cases. 

Figures  63  and  64  show  the  results  of  the  plane-wall  analysis.  It 
is  clear  from  a  comparison  of  Figs.  63  and  64  with  Figs.  57  and  58  that 
the  absence  of  the  body  gives  rise  to  significant  differences  in  the 
evolution  of  the  vortices.  Also,  the  dipping  of  the  tip  vortex  at  Z/bQ  =  14 
is  not  possible  without  the  streamlined  body.  Clearly,  the  half-body 
configuration  models  much  more  closely  the  whole-body  behavior  than  the 
plane-wall  model. 

The  foregoing  analysis  is  limited  by  such  assumptions  as  inviscid 
line  vortices,  instantaneous  vortex  roll-up,  no  boundary  layer  along  the 
body,  etc.  Nevertheless,  the  hydrodynamic  effects  of  the  body  are  strong 
enough  to  discourage  the  use  of  the  plane-wall  model  in  the  experiments. 
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Figure  62.  Real  and  Image  Vortices  for  the 
Plane  Wall  Analysis. 
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MIRROR  IMRGE  RNRLYSIS 
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Figure  63.  Positions  of  Tip  and  Root  Vortices  For  the  Plane  Wall  Analysis 
(Near  the  Lifting  Surface). 


TIP  VORTEX  POSITION 
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The  prediction  of  the  spiralling  of  the  tip  and  root  vortices  from 
the  analysis  of  the  whole  or  half  body  required  the  determination  of  the 
specific  conditions  under  which  the  two  said  vortices  cannot  spiral 
about  each  other.  Such  an  analysis  was  carried  out  only  for  the  simple 
case  shown  in  Fig.  65  where  the  body  effects  are  ignored.  The  use  of 
the  complex  function  for  the  two  real  vortices  and  their  images  together 
with  the  condition  of  no  spiralling  yielded 


a 


(E-7) 


where  a  3  jl^/Tgl  and  0  3  x^/x^.  Clearly,  there  exists  a  region  delineated 
by  certain  values  of  a  and  6  in  which  the  vortices  do  not  rotate  about 
each  other.  The  conditions  under  which  the  vortices  do  rotate  about  each 
other  may  be  exploited  for  more  rapid  demise  of  the  vortices.  However, 
it  remains  to  be  seen  as  to  how  and  under  which  conditions  the  root 
vortices  are  generated.  This  information  could  only  come  from  experiments 
and  is  very  much  related  to  the  lift  carrying  capacity  of  the  lifting 
surfaces  and  of  the  body  to  which  they  are  attached. 
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Figure  65,  The  Analysis  of  No-Rotation  Condition  of  the  Tip  and 
Root  Vortices. 
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